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ABSTRACT 


The depth of the marine atmospheric boundary layer (MABL) is an important 
parameter for both scientific and operational meteorological applications. The depth of 
the marine boundary layer has a significant influence on the atmospheric dynamics in the 
coastal zone. Knowledge of the depth of stratocumulus-topped boundary layers (STBLs) 
will enable coastal operations to more accurately anticipate weather, and electromagnetic 
propagation conditions. This study develops a satellite remote sensing technique for 
determining the height of MABLs topped with stratocumulus clouds. 

Validation of the technique using coastal rawinsonde dataset from the Monterey 
Area Ship Track (MAST) experiment revealed that an assumption of 41% cloud with a 
moist lapse rate equal to -7.0°C/km had the best overall fit to the data. However, for 
shallow boundary layers with depths below 400m the most accurate assumption was 75% 
cloud with a moist lapse rate equal to —-6.5°C/km. The application of this technique to 
sounding data returned an overall BL depth accuracy of 50m while the satellite 
application returned an overall accuracy of 65m. A sensitivity analysis of both surface 
and cloud-top temperature revealed that a 1/2°C change in either temperature resulted in 


an error of 60-70m in boundary layer depth. 
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I. INTRODUCTION 


The depth of the marine atmospheric boundary layer (MABL) is an important 
parameter for both scientific and operational meteorological applications. The depth of 
the marine boundary layer has a significant influence on the atmospheric dynamics in the 
coastal zone. Since the MABL is often topped by stratocumulus clouds, knowledge of 
the depth of stratocumulus-topped boundary layers (STBLs) will enable coastal 
operations to more accurately anticipate weather, and electromagnetic propagation 
conditions. A satellite remote sensing technique for determining the height of MABLs 


topped with stratocumulus clouds is developed in this thesis. 


a APPLICATION OF BOUNDARY LAYER DEPTH DATA 


The value of knowing boundary layer height can be seen when considering 
propagation conditions of surface based radar platforms. Even without complete 
knowledge of ducting profile shapes in the region knowledge of boundary layer height is 
important using the assumption that the boundary layer is well mixed. Boundary layer 
(BL) height is directly related to the base of the trapping layer associated with 
atmospheric ducts that can enhance the propagation of electromagnetic (EM) energy. The 
base of the trapping layer is referred to as the “optimum-coupling layer”. When 
conditions are such that the atmosphere acts as a wave-guide or duct, the propagation of 
EM waves is confined to a narrow region. Operational sensors that detect radar waves 
and other EM energy can operate more effectively when favorable ducting conditions are 


present. 





Surface vessels can experience enhanced EM detection conditions in the presence 
of surface-based ducts. The “optimum-coupling layer”, which is equivalent to the height 
of the well-mixed portion of the MABL, 1s related to the type of duct present. Figure 1-1 
illustrates the height difference for a surface-based duct (Z,) and for an elevated duct (Z;). 
The surface-based duct normally has a lower optimum-coupling layer (Davidson, 2000). 
The contrast between the two types of ducts and their impact on EM propagation 
conditions can be seen in Figure 1-2 which shows greater detection ranges for surface- 


based ducts vice elevated ducts (Rogers, L.T., 1999). 


Different Types of Ducts 





S00) eee: =. - Sa ee eee. - 
om 
= 
~ 200 

| 

Sb levated Duct Surface-based 
= Duct 
an 

100 

7a ae ee ee Evaporation Duct 


Figure 1-1: Examples of elevated and surface-based ducting 
conditions. 
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Figure 1-2: Radar image of elevated and surface-based ducting conditions. 


While boundary layer depth information is of critical importance for determining 
the type and effect of ducting conditions, the moisture content above the BL must 
also be know in order to determine if a trapping layer occurs, 1.e., to determine the 
shape of the refractivity (M) profile. Due to the presence of strong subsidence the 
air above the BL is generally very dry. Figure 1-3 shows a sounding from the 
Arabian Gulf, which clearly shows the dry nature of the atmosphere above the 
boundary layer in the presence of subsidence (Davidson, 2000). Ducting conditions 
are very common especially over eastern ocean basins where strong subsidence 
from a semi-permanent high is common. Figure 1-4 shows the importance of the 


knowledge and application of boundary layer depth to the determination of ducting 
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conditions can be seen when considering the prevalence of ducting conditions 


around the world. (Ortenburger et al, 1985) 
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Figure 1-3: Example of the nature of moisture profiles 


in and above the BL ducting in a region of strong 
subsidence. 


Rosenthal et al, [1997] showed that it is possible to obtain a statistical correlation 
of duct height (BL height) with cloud-top temperature (Ter) off the coast of 
California to within approximately 133m of accuracy 80% of the time (Rosenthal, 


J.S., 2000). The correlation of Ter to duct height implies that a remote sensing 


of 7 





technique for determining BL depth under strong subsidence inversions would 


provide accurate, real-time estimation of ducting conditions. 
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Figure 1-4: World wide frequency of ducting occurrence. (from 
Ortenburger, L.N., et al, 1985) 


B. METHODS OF CLOUD HEIGHT DETERMINATION 


Direct measurements of boundary layer height are often difficult and impractical to 
make and are therefore relatively rare. This has led to the development of several methods 
by which the height of the MABL may be directly and indirectly determined in cloudy and 
in cloud free atmospheres. In the past techniques for determining the height of the MABL 
were limited to land and ship launched radiosondes. However, several techniques have 


been developed to accomplish the same goal in a more efficient manner. Some examples 


of these techniques include ground-based SOund Detection And Ranging (SODAR) 
(Gaynor and Mandics, 1978) and LIght Detection and Ranging (LIDAR) (Eloranta et al., 
1975). These techniques require multiple stations and even in optimal circumstances they 
fail to provide a complete view of the marine boundary layer due to their limited spatial 
extent. 

The development of satellite-derived retrieval methods provides an expanded 
observation area and therefore a more complete picture of boundary layer depth. Many 
satellite techniques exist for measuring boundary layer properties including multispectral 
techniques for cloud free boundary layers (Kren, 1987). This study will focus on cloudy 
boundary layers since, as Kuze and Chance [1994] have shown, clouds cover more than 
50% of the Earth’s surface. 

There are several types of cloud height retrieval methods that utilize satellite data. 
These techniques can be broken up into two general categories. The first group extracts 
cloud top height in a direct way from satellite measurements including stereoscopic 
analysis and cloud shadow measurements. Stereoscopic analysis uses offset image pairs 
that are correlated to give a horizontal shift that 1s turned into a cloud top height. 
However, this technique is only applicable when near-coincident satellite image pairs are 
available for comparison and is therefore restricted to geostationary imagery (Shenk et al., 
1975). Cloud shadow measurements use the given solar and satellite geometry to calculate 
cloud height from the length and direction of observed cloud shadows. This technique 


applies only to clouds with limited vertical (Hasler et al., 1991). 


The second category uses the physical characteristics of the atmosphere to model 
cloud height. Some of these techniques include the use of atmospheric absorption bands 
(CO, or O), and bispectral estimations. Atmospheric absorption techniques use model 
optical path length and therefore height from within an absorbing region of the 
electromagnetic spectrum. Typical meteorological satellites such as the Advanced Very 
High Resolution Radiometer (AVHRR) instrument on the NOAA polar orbiters employs 
neither the correct wavelength regions nor the narrow bandwidths needed to estimate 
cloud top height. Bispectral methods on the other hand derive cloud-top height by 
comparing cloud-top brightness temperature (T3) with satellite derived vertical soundings 
such as from the Tiros Operational Vertical Sounder (TOVS). Unfortunately, satellite 
soundings lack the vertical resolution that is necessary for accurate estimations of STBL 
cloud height. Another option is to compare Tg with sounding output from numerical 
weather model temperature fields. These types of fields nominally have temperature errors 
of 1-2°C and this type of application generally yields cloud height errors of 300-500 
meters (Simpson, et al, 2000). 

The use of satellite data for broad spatial measurements of the depth of the STBL 
has the potential to significantly improve the understanding of large-scale STBL properties 
and dynamics. This study will present and test a technique for estimating the depth of the 
STBL by using a vertical thermodynamic model that calculates boundary layer depth from 
independent measurements of sea surface temperature (Ts) and satellite derived cloud-top 


brightness temperature. 


In Chapter II of this study the theoretical background behind the structure and 
dynamics of the STBL as well as the theory behind satellite derived cloud-top brightness 
temperatures will be presented. Chapter III will present the datasets utilized in this study 
as well as the specifics of the retrieval technique. Chapter IV will present the pertinent 
assumptions needed to accurately estimate the depth of the STBL and the results of the 
technique applied to atmospheric soundings and satellite derived Tg. Chapter V covers the 


conclusions drawn from this thesis research and recommendations for further study. 


Il. BACKGROUND 


With the proper assumptions, it is possible to indirectly determine the depth of the 
STBL using sea surface temperature and satellite derived cloud-top brightness 
temperature. The technique requires bnghtness temperatures from AVHRR channel 4 
(11pm wavelength), and an assumption of vertical cloud fraction in order to determine 
the height of transition from the dry to moist (pseudo-adiabatic) lapse rate within the 
STBL. This chapter covers the structure and dynamics of cloudy marine boundary layers, 
and the basic radiative transfer process that determines the measurement of cloud-top 


brightness temperature. 


A. LAPSE RATES IN THE STBL 


The Earth’s atmosphere cools with increasing height due to the fact that it is 
heated primarily from below by longwave radiation emitted from the Earth. The rate of 
cooling is known as the lapse rate of the atmosphere. The lapse rate of the atmosphere is 
associated with the mathematical slope of the line of temperature vs. altitude. By 
convention a decrease of temperature for increasing altitude is denoted as a positive lapse 
rate. 

The technique described here assumes adiabatic lapse rates within the STBL. 
However, the atmosphere near the Earth’s surface is adiabatic only under certain 
conditions. These conditions include heavy cloud cover, which limits radiative heating of 


the surface, and turbulence sufficient to mix the air and smooth out temperature gradients. 


In other words, the atmosphere may not be adiabatic due to insufficient mixing and/or 
solar heating of the surface. 

These two processes produce highly variable temperature behavior in the lowest 
layers of the atmosphere. Solar heating of the surface is most pronounced over land- 
covered areas. Over the ocean the high heat capacity of water prevents significant diumal 
warming and cooling of the ocean surface even during cloud free days. Therefore, the 
STBL can be said to meet these two critena most of the time, and an adiabatic lapse rate 
is therefore a good assumption. 


he Dry Adiabatic Lapse Rate 


The dry adiabatic lapse rate is the rate at which air temperature decreases with 
increasing altitude for an air parcel that is unsaturated. The decrease in temperature is 
due to the adiabatic expansion (occurring without loss or gain of heat) of a parcel of air. 
An expanding parcel of air will cool as the volume increases due to conservation of 
energy. The value for the dry adiabatic lapse rate is approximately -9.8°C/km. 


2. Moist Adiabatic Lapse Rate 


The moist adiabatic lapse rate is the loss of temperature with increasing altitude 
for an air parcel that is saturated. Ifa saturated air parcel is adiabatically expanded, the 
temperature will decrease, but this decrease causes the parcel to become increasingly 
Supersaturated. Saturation leads to condensation, and the latent heat released in this 
process warms the parcel. The warming from the latent heat release partially offsets the 
loss of temperature due to adiabatic expansion leading to a lapse rate that is smaller in 


magnitude than the dry adiabatic lapse rate. 
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It can be seen in Figure 2-1 that the pseudo-adiabatic lapse rate is a nonlinear 
function due to the addition of latent heat. However, this lapse rate is essentially constant 
below 1.5km. Since the STBL rarely exceeds 1.5km in height, a constant value is used 


for the pseudo-adiabatic lapse rate in this procedure. (Hsieh, 1987) 






y always stable 


Height (km) 


“es 
*. 
a 


unstable 


° 
“«, 
eo 
a 


-10 0 10 
Temperature (°C) 


Figure 2-1: Atmospheric lapse rates. 
(from Hsieh, 1987) 


B. STRUCTURE OF THE STBL 


For cloudless boundary layers the majority of the turbulence is produced at the 


surface by thermal convection and/or vertical wind shear. Surface turbulence leads to the 
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upward transport of momentum and moisture in the layer. In the STBL an additional 
source of turbulence is added by the presence of the cloud. Stratiform clouds at the top of 
the marine boundary layer result in dynamics that are significantly more complicated than 
for the cloud free boundary layer. The STBL is produced and maintained due to the 
balance of several processes. The nature of these processes will be discussed in the 
following sections. 


1 Cloud-Top Inversion 


STBLs form under inversion layers that are due to the descending, or subsiding air 
in the atmosphere associated with high-pressure areas. The inversion layer itself is a 
product of the adiabatic heating of the descending air and the layer just above the 
subsidence inversion 1s usually dry. Within an inversion layer there is a reversal of the 
normal lapse rate where temperature increases with altitude. The STBL is most prevalent 
beneath strong subsidence inversions over subtropical oceans. These areas include eastern 
subtropical ocean basins and mid to high latitude oceans where manine stratiform clouds 
are common. Specifically, the area in the North Pacific Ocean off the California coast is 
typically covered with marine stratocumulus clouds during the summer months when the 
region is dominated by a persistent sub-tropical high pressure system. 

The degree of turbulent mixing within the boundary layer determines the 
distribution of moisture, which influences the thickness of the cloud and can change the 
overall depth of the boundary layer. The following is a breakdown of the processes that 


produce the mixing required to maintain stratocumulus clouds. (Stull, 1988) 


Surface-based Free Convection: Free convection occurs when cold air is 
advected over a warmer surface since the condition of cold air over warm water 
results in an unstable atmosphere. The maintenance of stratocumulus clouds and 
fog off the coast of southern California is strongly influenced by surface-based 
free convection. 

Shear-generated Mechanical Turbulence: Shear turbulence is associated with 
strong BL winds that can cause convection (mixing) that is sufficiently strong 
enough to mix the moisture from the ocean surface to the cloud layer. The 
efficient redistribution of surface properties caused by mixing is an important 
factor in the maintenance of marine stratocumulus clouds. 

Cloud-top Radiative Cooling: The radiative heat loss at the top of 
stratocumulus clouds is a significant factor in the maintenance of the STBL. The 
cooling at the tops of stratocumulus clouds creates pockets of cold air that sink 
and mix with lower levels and help maintain the condensation necessary to 
produce cloud droplets. Cloud-top cooling is a process that is important to the 
maintenance of stratocumulus clouds as long as there are not higher clouds that 
can reduce the cooling of the lower cloud deck by radiating energy down onto the 
cloud top. 

Cloud-base Radiative Heating: The process of radiative heating, by radiation 
emitted from the surface and absorbed at the cloud-base, plays an important role 


in the destabilization of the cloud layer while stabilizing the subcloud layer. 
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However, it 1s generally much weaker than cloud-top cooling due to the smaller 
temperature difference between layers. 

Cloud-layer Shear: Shear generated in the cloud-layer generally produces 
mechanical turbulence and mixing near the top of the cloud. The turbulence 
resulting from cloud-layer shear produces small scale entrainment mixing that can 
decouple the cloud unless there are other turbulent processes present to more 
evenly distribute moisture and temperature throughout the BL. 


zs Cloud Base or LCL 


The altitude at which saturation occurs inside an air parcel, as 1t mses, expands and 
cools, is called the lifting condensation level (LCL). The LCL corresponds to the base of 
the clouds in the STBL. The value of the LCL is not constant in the STBL and its 
variation is governed by the degree of mixing or turbulence within the boundary layer and 
by a destructive process known as decoupling. 

Decoupling 1s the process by which stratocumulus clouds become disconnected 
from the lower regions of the marine boundary layer. When the boundary layer is 
decoupled a small inversion develops below the cloud base, which inhibits mixing from 
the ocean surface to the cloud and creates a lower, moist layer beneath a drier, cloud 
mixed layer. Decoupling occurs as a result of two primary processes in the boundary layer 
(Tjernstrém, 1998). 

a. Decoupling by Solar Absorption (Diurnal Decoupling) 
The process of diurnal decoupling is caused by variations in mixing due to 


cloud-top radiative cooling. During both day and night the longwave radiation that is 
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emitted upward from the cloud top is generally higher than the incoming energy received 
from the dry above the cloud. The net loss of radiative energy causes an imbalance that 
cools the top of the cloud layer and leads to thermally induced instability since the top of 
the cloud is now cooler than the bottom. The instability is enhanced under daytime 
heating as shortwave solar radiation penetrates into the cloud which warms it from 
within. 

The effect of solar radiation can cause a secondary inversion to form at the 
cloud base if the surface flux is weak. The secondary inversion separates the cloud from 
the surface and a second turbulent layer forms in the STBL. The diumal fluctuation of 
cloud thickness due to this process plays an important role in the understanding of the 


dynamics of the STBL. 


b. Decoupling Due to Precipitation 


The second decoupling process is due to precipitation out the bottom of 
the cloud. As described above, vertical cooling rate of the atmosphere is modified by the 
release of latent heat as water vapor condenses. However, as drizzle falls out of the cloud 
some of it evaporates in the air below leading to a cooler, moister sub-cloud layer. The 
development of a sub-cloud layer causes a stable layer to form near the cloud base that 
prevents the boundary layer from being well mixed. 


3. The Surface (Sea and Air Temperature) 


The surface air temperature (Ta) and sea surface temperature (T's) are usually 


close in magnitude. If a gradient in temperature occurs, fluxes develop which produce 
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differences in T, and Ts such that they seldom precisely equal in value. There are several 

physical process which can affect the differences in Ta and Ts. These processes include: 
Convection: Convection is a process by which heat is transferred due to the 
physical movement of air. Convective processes help to minimize temperature 
fluctuations by vertically mixing the air near the surface. 
Conduction: The transfer of heat through matter by communication of kinetic 
energy from particle to particle is known as conduction. Air is made up of many 
dispersed molecules that inhibit the process of conduction making air an effective 
insulator. The process of conduction is negligible when compared to the other 
process in the atmosphere, which are much larger in magnitude. 
Radiation: Radiation is the transfer of heat via electro-magnetic energy which is 
transferred to the atmosphere through absorption by gaseous constituents such as 
water vapor and carbon dioxide. Radiation from the ocean has a very strong 
moderating effect on the air temperature over the ocean’s surface, which is due to 
the large heat capacity of the ocean. The high heat capacity of the ocean enables it 
to gain and lose energy without large changes in Ts. 
Advection: In many cases, advection is the most influential of the four processes 
in amplifying the sea/air temperature difference. When an air mass moves out of 
its generation area and moves over ocean waters with different properties, strong 
fluxes of temperature and moisture result. These fluxes act to reduce the 


difference between the sea surface and the surface air temperature. 
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Klein and Hartmann [1993] showed the annual, climatological cycle of Ts and Ta 
difference to be very small in the Californian stratus region. Figure 2-2 shows the 
difference between T, and Ts is typically between 0.5°C and 1°C. This was done when 
they compared the seasonal variation of stratus cloud amounts off the coast of California 
using shipboard data. It 1s therefore assumed that the difference between the air and sea 
surface temperature 1s nominally within 1°C at any given time. This small (1°C) 
temperature difference is primarily due to climatologically weak surface fluxes of 
sensible heat. For simplicity in this study, it is initially assumed here that the sea surface 


temperature 1s equal to the surface air temperature. 
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Figure 2-2: The annual cycle of Ts and T, difference. 
(Klein and Hartmann, 1993) 


oF SATELLITE DETERMINED CLOUD-TOP TEMPERATURE 


The technique presented in this thesis uses the AVHRR sensor to determine 
cloud-top temperature. AVHRR channel 4 measures radiance at 11um wavelength. 
Channel 4 is ideal for measuring surface or cloud-top radiance due to the lack of gaseous 
absorption in the atmosphere at this wavelength. (Kidder and Vonder Haar, 1995) 


1. Radiative Transfer to Determine Cloud-Top Temperature 


It is assumed that the cloud layer being sensed is thick and that the atmosphere 
above the cloud top 1s dry. From Schwartzchild’s equation the radiance at the top of the 


atmosphere is: 


dt (A, p) 
L(A, 9,0) =, (A, 0)B(A, Ts)t, (A) + [ BAT) aw 


where L, is the spectral radiance at the top of the atmosphere as a function of wavelength 
(XA) and direction (9,@). €,(A,8)B(A,T;)tg(A) is the emitted radiance from the surface as a 
function of wavelength and surface temperature. tg(A) is the total direct transmittance 
throughout the full extent of the atmosphere as a function of wavelength and tg(A,p) is the 
direct transmittance from pressure p to the top of the atmosphere (p=0). B(A,T(p)) is the 
Planck blackbody emittance as a function of wavelength, and temperature at a given 
pressure. The derivative of the direct transmittance (tg) with respect to pressure 1s a 
weighting function that determines the contribution to L; due to at each vertical position 


in the atmosphere. 


For the STBL [dtg(A,p)/dp]dp = 1 over a few meter thick layer thick near the 
cloud top. Therefore, nearly all of the radiance emitted from the cloud comes from the 
upper portion only. The surface term [¢,(A,T;) ta(A)] goes to zero since the cloud is 
opaque at these wavelengths and it absorbs all of the energy from the surface. The final 


result is: 
L. = B(A,Te,) (2-2) 


Here L, is measured by the sensor at A and the value of Ter can be readily determined 
from inversion of the Planck blackbody relation. 


Ze Errors in Satellite Determined Tcr7 of the STBL 


Several conditions can potentially interfere with the determination of Tc; from 
satellite Tp. 
a. High Clouds 
The presence of higher clouds (e.g., cirrus) above the STBL produces 
lower values of Tcr due to the lower emitted radiance of clouds that exist in the higher, 
colder levels of the atmosphere. The presence of cirrus clouds in particular can be a 
problem due to the fact that cirrus ones are often thin and difficult to detect without 
multispectral techniques. Therefore, the presence of cirrus clouds in the vicinity of the 
STBL can introduce significant errors in low level, satellite derived cloud-top 


temperatures. 
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b. Limb Darkening 

When a measurement is taken from a satellite at the periphery of its field 
of view the radiance that 1s detected must pass through a larger volume of atmosphere 
than radiance measured from nadir. The effect of greater atmospheric interaction reduces 
the amount of cloud-top radiance that reaches the satellite sensor. The AVHRR sensor 
has a scan angle from nadir of +/-55.3° (Kidder and Vonder Haar, 1995). Measurements 
made at an angle of 55.3° are not only made through an increased atmospheric volume 
but radiance emitted from the cloud-top has an angular dependence as well. The radiance 
emitted from a cloud decreases as the view angle increases. Therefore measurements of 


Tcr made at the edge of a satellite pass are often significantly reduced in radiance. 


D. THE VALIDITY OF A WELL-MIXED STBL 


Within the STBL the potential temperature and water vapor mixing ratio are 
homogeneous due to turbulent mixing within the layer. The eddies which produce this 
mixing are due to the presence of buoyancy and shear effects. The magnitude of 
buoyancy production within the STBL is small due to the high heat capacity of the ocean 
which prevents large diurnal temperature fluctuations. Nighttime radiative cooling of the 
sea surface prevents the establishment of a near surface stable layer. Therefore, the STBL 
is often neutrally stable with only slight diurnal variations. (Kren, 1987) 

The weakly unstable or neutral STBL is influenced by wind shear as well. The 
effect of wind shear on the STBL is to produce turbulent eddies which mix the STBL. 


Turbulence produced by both buoyancy and shear is generally sufficient to mix the STBL. 
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Some regions under synoptically induced subsidence inversions have a stable 
STBL where mixing is often suppressed. These areas have relatively cold sea surface 
temperatures that often comes from the coastal upwelling found along the west coast of 
continents. In these areas where the STBL is not well mixed, using a dry adiabatic lapse 
rate will add uncertainty to the final value of boundary layer depth. The accuracy of the 
technique used in this thesis under these conditions will be described in the next few 
chapters of this text. 

The technique presented in this thesis has certain limitations which constrain its 
applicability. The technique requires knowledge of the temperature of both the sea 
surface and the cloud-top. However, satellite techniques for obtaining sea-surface 
temperature require a cloud-free atmosphere, and techniques for obtaining cloud-top 
temperature requires either an atmosphere with little water vapor above the STBL, or 
knowledge of the water vapor content so its effects can be corrected. The technique 
presented here has potential application in regions of strong subsidence inversions where 
cloudy boundary layers persist. Sea surface temperature data can be provided by ship 
observations, recent satellite analyses of SST from clear sky periods or predictive ocean 


models. 
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Iii, DATA AND METHODOLOGY 


ae MONTEREY AREA SHIP TRACK (MAST) EXPERIMENT DATASET 


The primary data set for this thesis comes from the Monterey Area Ship Track 
Experiment (MAST) (Durkee et al, 2000). MAST was a multi-agency effort designed to 
determine the physical nature of ship tracks via the collection of sounding, aircraft and 
satellite data off the coast of Central California. Figure 3-1 depicts the MAST operating 
area. This thesis will focus on AVHRR satellite data and radiosonde data collected by the 


R/V Glorita during the MAST experiment. 





Figure 3-1: MAST Experiment 1994 Operations Area. 
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1. Radiosonde Data 


The R/V Glorita launched an average of six radiosondes a day at intervals of 
roughly four hours for a total of 94 radiosonde soundings (designated GLO1-GL94). 
Parameters measured with the radiosondes were pressure, temperature and relative 
humidity. These measured quantities where then used to derive values for height, 
dewpoint and mixing ratio. The techniques used to derive these quantities are included in 
Appendix B. 


2: Surface Observation Data 


Personnel onboard the R/V Glorita made hourly measurements of sea surface 
temperature (Ts), air temperature (T,), relative humidity, as well as wind speed and 
direction. However, gaps in Ts measurements did occur when the temperature sensor was 
not in the water. Additionally these measurements are amplified by general comments on 
sky conditions made during radiosonde launch times. 


3. Satellite (AVHRR) Data 


Imagery from AVHRR Channels 1-5 of the NOAA-9, 10, 11, and 12 satellites was 
recorded during the experiment. Channel 4 (11um) brightness temperatures from this 
imagery were used to determine the cloud top temperature at the location of a given 


sounding from the R/V Glorita. 


B. ANALYSIS PROCEDURE 


In order to determine the depth of the STBL a vertical, thermodynamic model was 


constructed. The model assumes that unsaturated air follows the dry adiabatic lapse rate 
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(Ig) until the LCL is reached. From the LCL to the cloud top the air is saturated and 
follows a pseudo or moist adiabatic lapse rate ([,). With a given assumption of cloud 
thickness these two lapse rates can be combined in order to give a more complete 
thermodynamic picture of the STBL. 

The combination of Ty and I, determines a cloud-top height (Z) for a given 
surface temperature (Ts) and cloud-top temperature(Tcr). The constant lapse rate 
equation I = (Tcr-Ts)/Z is simply rearranged algebraically to yield the formula for height 
Z=(Tcr-Ts)/I. Given a formula for height based on temperature and lapse rate only, it is 
now possible to compute heights based on the moist and dry adiabatic lapse rates for 
given values of Tcr and Ts. 

An assumption of vertical cloud fraction must be made in order to give a proper 
physical parameterization of how much of the boundary layer follows I, and how much 
follows Ig. Since the actual height of the boundary layer is unknown an initial value, or 
first guess, for the boundary layer height must be used in order to determine the LCL or 
cloud base height. The LCL in this technique 1s derived from an assumed fraction of the 
total, dry adiabatic height (Zy). 

Figure 3-2 shows an example where it is assumed that the upper one-third of the 
boundary layer is cloud and the remaining two-thirds of the boundary layer is cloud-free 
air. The LCL is determined by multiplying the total height of the cloud-free boundary 
layer (obtained via the dry adiabatic lapse rate) by two-thirds. The remainder of the 


boundary layer height is calculated from the LCL by using the pseudo-adiabatic lapse rate 
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until the value of Tcr is reached. However, by virtue of the fact that the pseudo-adiabatic 
lapse rate is steeper than the dry lapse rate this process produces a slightly larger fraction 
of cloud than is initially assumed. Therefore, when an assumption of 1/3" cloud is made 
based on the height of a totally dry BL, the increase in height from the moist-adiabatic 


lapse rate produces a final BL distribution of 41% cloud and 59% cloud-free atmosphere. 
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Figure 3-2: The application of the assumption of a cloud thickness of 
41% and how it 1s derived from the dry adiabatic height Zd. 


In order to calculate the height of the lower, cloud-free portion of the boundary 
layer the height of the entire boundary aval (i.e., from Ts to Tcr) must be computed 
along the dry adiabat with the relation Zg = (Tct-Ts)/ Ig where Z, 1s the height of the 
boundary layer if it were entirely cloud free. The height of the cloud free portion of the 


boundary layer is then simply the height of the cloud-free BL multiplied by a slight 
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overestimate of the initial fraction of cloud-free STBL or Zcr = Zg*(Fcr). Where Zcr 
denotes the initial estimate of the height of the cloud-free STBL and Fer is the initial 
estimate of the cloud-free fraction of the STBL. 

Once the initial height of the cloud free portion of the boundary layer (which is 
actually the height of the cloud base) is determined, the corresponding temperature at that 
height must be known in order to calculate the height of the saturated portion of the 
boundary layer. The temperature at the cloud base or LCL is computed by the relation 
Tcr = Zcr*lg + Ts which then allows for the same fundamental lapse rate relation to be 
used in the upper two thirds of the boundary layer. This final height is combined with the 
height of the below cloud BL to yield the total height of the STBL via the following 
relation: Zp = Zcr + (Tct-Tcer)/ Im. 

In order to determine the accuracy of the technique, values of T, and Tcr were 
taken directly from the Glorita soundings and used to determine the height of the STBL. 
To determine the accuracy of the technique using satellite data, the inversion temperature 
from the sounding is replaced by the cloud-top brightness temperature from the AVHRR 
channel 4 imagery. A total of 46 soundings were used from the MAST dataset and 48 
soundings were dismissed outright from the MAST dataset due to insufficient MABL 
clouds. See Table 3 in Appendix A for more detail on which soundings were dismissed 


from the dataset. 
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IV. RESULTS 


A. ASSUMPTIONS 

Two primary assumptions must be made in order to model the general state of the 
STBL. The amount of cloud must be specified and the lapse rates of the atmosphere in 
the cloud and below cloud must be given. In the STBL it is assumed that all lapse rates 
are of constant value and near adiabatic in nature. The below cloud atmosphere is 
assumed to be exactly dry adiabatic with a lapse rate of gy = -9.8°C/km and within cloud 
the value of the pesudo-adiabatic lapse rate is initially given the standard atmospheric 
value of [y= -6.5°C/km. 

It was initially assumed that the most common ratio of saturated to unsaturated air 
in the STBL was 2/3™ cloud to 1/3™ unsaturated (below cloud) air. However, results 
have shown that a more accurate depiction of the STBL requires an assumption of 75% 
cloud for shallow boundary layers (below 400m deep) and a value of 41% cloud for 
boundary layers deeper than 400m. Above 400m the value for the pseudo-adiabatic lapse 
rate was modified to a value of [m= -7.0°C/km to partially compensate for its non-linear 
nature as altitude increases. The reasons these assumptions were chosen will be 


discussed below. 


B. SOUNDING APPLICATION 


In order to validate the physical reasoning of the technique described in Chapter 


Ill, the method was applied to the soundings from the R/V Glorita. The method was run 
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with several different values for I,, and cloud thickness in order to determine the best set 
of parameter values. The results for each set of assumptions were organized by the 
difference between the actual BL height and the calculated height of the STBL (denoted 
as AH). Throughout the venfication process the assumptions made about the STBL were 
modified in order to minimize the value of AH 


2 Actual Depth Method 


The distribution of AH with respect to the actual BL depth was plotted in Figure 
4-1, which shows how the two sets of assumptions (75% and 41% cloud) begin to diverge 
in the vicinity of 400 meters. Based on these criteria it was determined from the sounding 
results that there was a discontinuity in the accuracy of a single set of assumptions at 
around the 400m point. Boundary layer cases with depth less than 400m are more 
accurately modeled by an assumption of 75% cloud with a lapse rate of [,, = -6.5°C/km. 
These values were the only set of assumptions that returned a value for AH < 100m in all 
cases below 400m. For boundary layers deeper than 400m it was determined that an 
assumption of 41% cloud combined with a lapse rate of IT, = -7.0°C/km returned the 
most accurate height. These assumptions for soundings above 400m produced a height 
differential (AH) that was consistently below 100m for the majority of the soundings. 

Figure 4-2 shows how the two sets of assumptions were combined into a two- 
layer model in order to give a more accurate picture of the total STBL. The height 


produced by this two-layer procedure were plotted against the actual BL height with the 
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Figure 4-1: Scatter plot of the two primary 
sets of assumptions applied all soundings. 


The results of this regression produced a line with a slope of 1.0283 where a slope 
of 1.0 is a perfect correlation. The line was offset from the origin by a value of -41m 
while the overall error of the estimate of this technique is 56m. Table | in Appendix A 
shows how several of the different assumption sets affect the determination of boundary 


layer height. 
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Figure 4-2: Linear regression of actual vs. calculated 
boundary layer height using sounding data. The solid line 


is the regression line and N is the number of data points. 


De First Guess Method 


While this process produces favorable results it requires previous knowledge of 
the actual depth of the STBL. Since the actual depth of the STBL is the goal of this thesis 
it cannot be used experimentally to determine the dividing line between the two sets of 
assumptions. Therefore, a “first guess” must be made in order to determine which set of 
assumptions to use. 

It can be seen in Table 1 that for a single set of assumptions the most accurate 
parameters over the entire STBL are a 41% cloud layer and pseudo-adiabatic lapse rate of 
Tm =-7.0°C/km. These assumptions were applied to the entire STBL and once an initial 


value of BL depth was obtained, the “shallow” set of assumptions were applied to the 
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soundings below the depth of 400m. Figure 4-3 shows the results from this two-step 
process where soundings were plotted against the actual BL depth and a linear regression 
was performed. The regression was found to be even more accurate than the original 
application of the assumptions with a slope of 1.0080, an offset of -24m and an overall 
error estimate of 50m. 
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Figure 4-3: Linear regression of actual vs. calculated 
boundary layer height using sounding data. The solid line is 
the regression line and N is the number of data points. 
* _ Indicates the composite BL height was determined by 
the two-step process of the application of BL assumptions. 
The soundings from the R/V Glorita indicate that shallow BLs, (Zp. < 400m) 
have a larger fraction of the BL that is saturated (cloudy) air than for deeper BLs. 


However, the height value of 400m is by no means an absolute dividing line between BLs 


with differing percentages of cloud. The MAST dataset contains several cases where 
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soundings with BL depths above 400m contain less cloud than 75%. When the 
assumption of 41% cloud is applied to these cases a significant underestimate is returned. 
This occurs due to the fact that when the assumption of 41% cloud is applied to these 
cases, the predominance of the dry lapse rate produces a significant underestimate in BL 
depth. When these soundings are reevaluated with the shallow assumption set (75% 
cloud), a more accurate determination of the depth of the STBL results due to the steeper 
pseudo-adiabatic lapse rate. Therefore, the process of making a first guess at the depth of 
the BL transforms the dividing point for the two sets of assumptions into a broader range 
of depths. For this dataset the effect was to extend the depth from 400m up to 460m in 
cases where the BL depth was severely underestimated by the assumption of 41% cloud. 


3. Sounding Case Studies 


Table 3 lists all R/V Glorita soundings along with position, launch time and 
general comments on their analysis. A total of 46 soundings were studied in the sounding 
application part of this thesis in order to validate the various assumptions made about the 
properties of the STBL. The following is a breakdown of several soundings that are 
representative of three primary depth categones. These cases are presented in detail in 
order to demonstrate how the assumptions described above apply to the various boundary 


layer depths encountered during this study. 


a. GL16 


Figure 4-4 shows a plot of sounding GL16 that was launched at 1153 UTC 


on 11 June 94 at 36.15°N-122.57°W. The surface pressure for this sounding is 1015.0 
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mb, relative humidity is 98.8%, and winds are from the NNW at 0.3m/s and increasing in 
strength with height. Ts; was not recorded for this case and surface air temperature is 


13.2°C. The depth of the boundary layer in this case is 330 meters with a cloud fraction 
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Figure 4-4: Sounding GL16 with lapse rate assumptions of 
75% cloud and I, = -6.5°C/km applied. 


When the procedure is run on the sounding GLI16 the first guess 
assumption (41% cloud) returns a BL height of 254m which corresponds to a AH of — 
76m or an underestimate of 76m. This falls in the <400m category and the technique 
requires a shift in the assumption set to the 75% cloud assumption. This combination of 
lapse rates gives a much more favorable result where the calculated height is now 300m 


with a AH of only —30m. 
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b. GL67 

Figure 4-5 shows a plot of Sounding GL67 that was launched at 1148 
UTC on 22 June 94 at 35.78°N-123.18°W. The surface pressure for this sounding is 
1019.8 mb, relative humidity is 77.6%, and winds are from the NW at 9.6m/s. Surface air 
temperature is 14.1°C, but Ts was not recorded for this case. The depth of the boundary 
layer in this case 1s 504 meters. 

When the procedure is applied to sounding GL67 the first guess 
assumption of 41% cloud returns a height value that is well above the 400m height and 
therefore these values will be retained in dealing with this sounding. The 41% cloud 
assumption returned a calculated height of 507m with a AH of +3m. 

It 1s important to note that when this sounding 1s compared to a satellite 
image that is described later it can be seen that it passes through an area of reduced 
cloudiness that is immediately adjacent to a broad area of very even marine stratocumulus 
clouds. Due to the uniform nature of the surrounding cloudiness it 1s assumed that the 
boundary layer properties associated with the surrounding clouds is represented in the 


sounding GL67 even though the sounding does not pass directly through these clouds. 
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Figure 4-5: Sounding GL67 with lapse rate assumptions of 
41% cloud and I’, = -7.0°C/km applied. 


C. GL68 


Figure 4-6 shows a plot of Sounding GL68 that was launched at 1453 UTC 
on 22 June 94 at 35.70°N-123.46°W. The surface pressure for this sounding is 1021.1 
mb, relative humidity is 84.4%, and winds are from the NW at 9.0m/s and decreasing in 
strength with height. The surface air/sea temperature difference is T,-Ts = -1.1°C, and 
the surface air temperature is 14.0°C. The depth of the boundary layer in this case is 


671m with a cloud fraction of approximately 37%. 
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Figure 4-6: Sounding GL68 with lapse rate assumptions 
of 41% cloud and I, = -7.0°C/km applied. 


When the lapse rate procedure is applied to this sounding the first guess 
assumption is utilized for this case due to the returned boundary layer depth of 665m, 
which corresponds to a AH of —-6m. As a test case the actual Ts value (Ts = 15.1°C) was 
used in this method and the result was an increase in the calculated depth of the boundary 
layer. This initial temperature modified the result such that an assumption of 41% cloud 
returns a calculated height of 789m with a AH of +118m. This indicates that a 1°C 
change in surface temperature modifies the height differential value by approximately 
120m. The effect of surface and cloud-top temperature variations on the entire sounding 


dataset will be discussed in section C of this chapter. 
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d. GL17 


Figure 4-7 shows a plot of Sounding GL17 that was launched at 1448 UTC 
on 11 June 94 at 36.23°N-122.48°W. The surface pressure for this sounding is 1015.5 


mb, relative humidity is 99.3%, and winds are from the West at 1.0m/s. 
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Figure 4-7: Sounding GL17 with lapse rate assumptions of 
75% cloud and I, = -6.5°C/km applied. 


The surface air temperature for this sounding is 12.2°C and the value of Ts was not 
recorded for this sounding. The depth of the boundary layer in this case is 340m with a 
cloud fraction of approximately 88%. 

Sounding GLI17 is relatively shallow and when the first guess set of 


assumptions are applied it returns a height of 242m and therefore the 75% cloud 
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assumption 1s applied to this sounding. The assumption of 75% cloud gives a calculated 


BL height of 286m which corresponds to an underestimate of AH =—54m. 


e. GL29 


Figure 4-8 shows a plot of Sounding GL29 that was launched at 2347 
UTC on 12 June 94 at 37.05°N-123.25°W. The surface pressure for this sounding is 
1016.8 mb, relative humidity is 91.9%, and winds are from the NW at 7.0m/s and 
increasing in strength with height. The surface air temperature for this sounding is 14.0°C 
and the air/sea temperature difference 1s Ta-Ts = -0.6°C. The depth of the boundary layer 
in this case 1s 377 meters with a cloud fraction of approximately 23.5%. 

This sounding initially meets the conditions necessary to be reevaluated 
with the assumptions for BLs below 400m. The first guess routine retums a BL depth of 
277m with a corresponding height differential of AH =-100m. When the assumption of 
75% cloud is applied to this sounding is returns a boundary layer height of 327m which 


corresponds to a height differential of AH = -50m. 
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Figure 4-8: Sounding GL29 with lapse rate assumptions of 75% cloud 
and I, = -6.5°C/km applied. 


C: TEMPERATURE VARIATIONS 


It was shown in the previous section that variations in surface temperature can 
have a significant effect on the accuracy of the technique presented in this thesis. Figure 
4-9 shows an individual example where sounding GL29 was modified such that the lapse 


rate is started at the actual sea surface temperature (Ts = 14.6°C). 
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Figure 4-9: Sounding GL29 corrected to start at the actual SST of 14.6°C. 
The steep, light blue lines represent the two different values for the 
pseudo-adiabatic lapse rate. 

A sensitivity test was made for the entire sounding dataset where both surface and 
cloud-top temperature were modified. The results of these modifications are compared to 
the 41% cloud assumption as well as a “no cloud” and an “all cloud” state. Figures 4-10 
thru 4-13 demonstrate how an offset in surface and cloud-top temperature (Ts, Tr) affects 
the accuracy of the model. Figures 4-10 and 4-12 indicate that a '2°C decrease in Ts or a 


/2°C increase in Tez produces a resulting accuracy that is equivalent to that of a model that 


relies solely on a dry adiabatic lapse rate for the entire BL. That is, these temperature 
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changes approximate that of a model that utilizes a dry adiabatic lapse rate over the entire 
depth of the BL. Figures 4-11 and 4-12 show how a 1°C increase in Ts or a 1°C decrease 
in Ter produces a result that is equivalent to that of a model that relies solely on a dry 
adiabatic lapse rate for the entire BL, 1.e., a boundary layer completely filled with cloud. 
The blue (41% cloud) line as with all of the lines represents a single set of 
assumptions applied to the entire dataset. The slight increase in error or tail that lies away 
from and below the main peak represents the shallow cases where Zp, <400m. These 
figures give a sense of the impact of inaccuracies in both Ts and Tc; where even relatively 
small errors in temperature can have a significant impact on the accuracy of the technique. 
These results imply that if errors Tc; or Ts are greater than 1°C then there is essentially no 


advantage to breaking up the lapse rates when determining BL depth. 
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Figure 4-10: The effect of a decrease in surface temperature compared 
to a cloud-free boundary layer. 
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Figure 4-11: The effect of an increase in surface temperature 
compared to an all-cloud boundary layer. 
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Figure 4-12: The effect of an increase in cloud-top temperature 
compared to a cloud-free-boundary layer. 
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Figure 4-13: The effect of a decrease in cloud-top temperature 
compared to cloud-free boundary layer. 
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D. SATELLITE APPLICATION 


The sounding profiles from the Glorita were paired with images from NOAA- 
9,10,11 and 12 that were closest in time to the soundings. Due to the relatively 
homogeneous nature of environment that produces marine stratocumulus clouds a 
maximum time difference of three hours was chosen as the criteria to identify 
sounding/image pairs. 

Channel 4 of each satellite image was analyzed in the region that the sounding 
was taken in order to determine the cloud top temperature. However, the presence of 
cirrus clouds in the image caused some cases to be dismissed. Cirrus clouds, being much 
colder than MABL clouds, contaminated the bnghtness temperature and produce 
erroneous results that significantly overestimate the depth of the STBL. Therefore, this 
method is not valid for images with significant cirrus contamination. For detailed 
information on which soundings were utilized for this procedure and which were 
dismissed see Table 4 in Appendix A. 

The technique applied to the sounding dataset was applied to satellite imagery in 
order to evaluate the accuracy of the application of Tg measurements that are utilized as 
cloud-top temperature. The value of Ts in this part of the study was once again obtained 
from the sounding in order to better understand the accuracy of substituting Tg for Tcr. 
The returned values of BL height were plotted against the actual BL height with the actual 
BL height as the independent variable. 

Figures 4-14 and 4-15 show a linear regression was performed on the data points 


in order to determine the accuracy of the process for both the direct method and the first 
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guess method respectively. The result from the first guess method was again found to be 
even more accurate than the direct (single assumption set) application. <A _ linear 
regression of the direct method returned a slope of 1.3074, an offset of -153m and an 
overall error estimate of 76m. The regression of the first guess method returned a line 
with a slope of 1.2687, an offset of -137m and an overall error of 65m. The details of 
this application and how it compared to the ground truth of the soundings can be seen 


Table 2 in Appendix A. 
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Figure 4-14: Linear regression of actual vs. 
calculated boundary layer height using satellite 
derived T, data. The solid line is the regression line. 
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Figure 4-15: Linear regression of actual vs. 
calculated boundary layer height using satellite T, 


data. The solid line is the regression line. 
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Ik. Favorable Satellite Case Studies 


The accuracy of using satellite Tgp as cloud-top temperature in the technique 


described above will be illustrated below for three sounding image pairs. 


a. GL16 


Figure 4-16 shows the closest satellite image associated with 
sounding GL16 that was taken at 2349 UTC on 27 June 94. The visual image in 


the area of this sounding shows a uniform sheet of stratocumulus cloud. 
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Figure 4-16: AVHRR Channei 2 image for GL16. 


Figure 4-17 shows the corresponding Tg taken from the IR mage has a value of 
11.4°C, which is very close in magnitude to the sounding cloud-top temperature 
(11.0°C). The application of this Ts to the lapse rate calculation of BL height 


yields a height of 251m (AH = -79m). 
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Figure 4-17: AVHRR channel 4 image for GLI1O6. 


For comparison the sounding based technique had a height differential of 


AH = -30m showing that, for this situation, a 49m change in accuracy resulted from 
a Ter error of +0.4°C. 
b. GL67 
Figure 4-18 shows the image that most closely corresponds to sounding 
GL67 was taken at 1148 UTC on 22Jun94. The corresponding Ts taken from the IR 
§ S 4B 


has a value of 10.0°C, which is very close in magnitude to the sounding cloud-top 


image 
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Figure 4-18: AVHRR Channel 4 zoomed image from for GL67 


The sounding in this case appears to go through a region of thin cloud the details of which 
can be clearly seen in Figure 4-18. The corresponding value of Ts measured in the IR 
image was obtained from the nearby thicker cloud which were assumed to be 
representative of the BL in the vicinity of sounding GL67. The application of this Ts to 
the lapse rate calculation of BL height vields a height of 477m (AH = -29m). For 
comparison the sounding based technique had a height differential of AH =+3m showing 


that. for sounding GL67, a 32m change in accuracy resulted from a Tcr error of +0.3°C. 


roa GL68 
The image that most closely corresponds to sounding GL68 was taken at 
1514 UTC on 22Jun94. The visual image (Figure 4-19) in the area of this sounding 


shows a region of even, well distributed stratocumulus clouds. 





Figure 4-19: AVHRR channel 2 image for GL68. 
The corresponding Tz taken from the IR image (Figure 4-20) has a value of 8.4°C. which 
is very close in magnitude to the sounding cloud-top temperature (8.2°C). The 
application of this Tg to the lapse rate calculation of BL height vields a height of 647m 
(AH = -24m). For comparison the sounding based technique had 2 height differential of 
AH = -6m showing that, for this situation. a 18m change in accuracy resulted from a Ter 


error of ~0.2°C. 
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Figure 4-20: AVHRR channel 4 image for GL68. 


ple Unfavorable Satellite Case Studies 


There are several situations where the accurate retrieval of Ter from satellite 
imagery is not possible. Errors in Tcr can occur when there are higher clouds (e.g., 
cirrus) above the STBL. or when there 1s a shallow BL where the clouds are too thin and 
permit radiance from the surface to contaminate the Tg readings. The latter case will be 
illustrated below as will an instance of “limb darkening” where Tg are obtamed from the 


far edges of the satellite path. 
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a. GLI7 
The closest satellite image associated with sounding GL17 was taken at 
1531 UTC on 11 June 94. This image can be seen on Figures 4-21 and 4-22. which show 


the visual (channel 2) and the infrared (channel 4) respectively. 





Figure 4-21: AVHRR channel 2 image for GL17. 


The brightness temperature derived from channel 4 for this mage was 10.8°C, which 1s 
slightly warmer than the cloud-top temperature measured on the sounding (10.1°C). This 
increase in Tg is possible due to the influence of continenta! aerosols that reduce droplet 
size and absorptance in the stratocumulus clouds leading to radiance from lower in the 


cloud contributing to Tg and producing higher values of Tg (Brenner. J.R.. 1994). 


34 


The sounding technique results that were described previously in this 
chapter returned a BL height of 286m or an underestimate of AH = -64m. The 
application of the warmer Tp from the satellite image produced an enhancement of this 
underestimate and returned a BL height of 192m ora AH = -157m. This indicates that 


there was a 93m change in accuracy due to a +0.7°C change in T¢r. 
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Figure 4-22: AVHRR channel 4 image for GL17. 
b. GL8&2 
Sounding GL82 was launched at 2349 UTC on 27 June 94 at 35.17°N- 
124.82°W. The surface pressure for this sounding is 1021.8mb, relative humidity ts 


83.8%, and winds are from the NNW at 9.2m/s. The surface air/sea temperature 
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difference is Ta-Ts = 0.0°C, and the surface air temperature is 16.1°C. The depth of the 
boundary laver in this case 1s 640m with a cloud fraction of approximatelv 34.3%. 
The closest satellite image associated with sounding GL82 was taken at 


2334 UTC on 27 June 94. Figure 4-23 shows the IR satellite image for GL82. The image 


shows that the corresponding sounding position was on the far left-hand side of the 
swath. This situation falls into the “limb darkening” category where radiance from the 


position of the sounding is reduced in intensity due to the extreme angle at which it 


travels to reach the satellite sensor. 





Figure 4-23: AVHRR channel 4 image for GL82. 


The channel 4 brightness temperature measured at the location of the 


| arta nie a e T , Se eles hich 3 : 7 tly ralder than ths 
unding on Figure 4-23 has a value of 9.1°C, which is significanily colder than tne 
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corresponding cloud-top temperature taken from sounding GL82 (Tc; = 10.2°C). This 
temperature differential results in a calculated BL depth of 808m which corresponds to a 
height differential of AH = +168m. 

3. Large Scale Satellite Application 

The next step in the logical progression of this technique is to apply it to the entire 
image vice a single point. This was done using a single Ts measurement from sounding 
GL62 to represent the entire Ts field along with the corresponding AVHRR image shown 
in Figures 2-24 and 4-25 that were taken at 1535 UTC on 21 June 94. The result of this 
process can be seen in Figure 4-26 which represents a topographic map of BL depth off 
the coast of Califorma. The image shows the depth of the boundary layer increasing away 
from the coast, which agrees with existing theones on boundary layer depth near 
coastlines (Stull, R.B., 1988). Depth (in meters) can be determined on the image by use 
of the color scale located on the upper left-hand side of the image. The image contained 
higher clouds (seen in the contrast between Figures 24 and 25) and cloud free areas which 


were filtered out and appear black on the image. 
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Figure 4-24: AVHRR channel 2 image for sounding GL62. 
| : 2 Jun 94 
1535 UTC 
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Figure 4-25: AVHRR channel 4 Enhanced IR image of GL62. 
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Figure 4-26: Topographic map of BL depth based on sounding GL62’s spel 
and AVHRR channel 4 imagery. The red cross indicates the location of 
sounding GL62. 
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ie CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 


The purpose of this thesis was to develop and test a technique for determining the 
depth of the STBL using boundary layer thermodynamics. The technique determines BL 
depth by utilizing satellite derived cloud-top brightness temperature and sea surface 
temperature observations. The method was initially refined with sounding data from the 
MAST experiment which were used to verify the technique. The technique was then 
applied to satellite brightness temperature data in place of the cloud-top temperature from 
the soundings. 

Analysis of the MAST dataset revealed that an assumption of 41% cloud (59% 
cloud-free) with a moist lapse rate of I’, = —7.0°C/km produced the best overall fit to the 
data. However, for shallow boundary layers the most accurate assumption set was 75% 
cloud with a lapse rate of I’, = -6.5°C/km. The shallow cases were found to apply to 
boundary layers with depths less than 400m. The value of 400m was chosen because it 
was found to be the center of a gradual shift in the accuracy of the two assumption sets. 

The sensitivity of the technique was analyzed by varying the values of both 
surface, and cloud-top temperatures. Analysis revealed that a '2°C change in either 
temperature resulted in an error in calculated BL depth of 60-70m. The magnitude of BL 
height error was found to vary linearly due to the linearity of the lapse rates used in the 
calculation of BL depth. The sensitivity of this technique to temperature inputs implies a 


strong need for accurate T; and T¢7 fields whenever this technique is used. 
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This technique is applicable to regions where the atmosphere above the BL is 
essentially cloud free as cirrus clouds significantly impair the accurate determination of Tg 
from the STBL. Furthermore, the technique presented here was designed for coupled BLs 
in a uniform cloud field. This technique will not yield accurate results if applied to 


decoupled BLs. 


B. RECOMMENDATIONS 


This study quantitatively demonstrates the sensitivity of BL depth determination to 
the accuracy of surface and cloud-top temperature data. Models that generate surface 
temperature data are generally of much lower resolution than the 1km accuracy of 
AVHRR Ts data. Future work in this area should ensure that the resolution of both Ts; and 
Ter fields are compatible. The application of Ts data to the procedure presented in this 
study will introduce errors due to differences between Ts and surface air temperature. The 
errors introduced into the technique due to this difference could be minimized if surface 
air temperature was obtained from a coupled model that could account for air/sea 
temperature differences. 

Another area of refinement is in the sensitivity of remotely sensed cloud-top 
temperature to cloud thickness. Measurements of Tg taken from thin clouds can be 
contaminated with surface radiance, which makes the value of Tg too high. Multi-channel 
remote sensing techniques can be employed to estimate cloud optical thickness and correct 


the measured Tp to account for contamination from radiance below the cloud-top. 
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The application of this technique to EM propagation and ducting conditions can be 
refined with the application of multi-spectral remote sensing techniques. Multi-spectral 
techniques could be employed to enhance the determination of the moisture content above 
the boundary layer which will help with the determination of refractivity conditions (i.e., 
ducting). 

Additional study should be done with this technique in order to better determine 
the nature of shallow BLs (<400m) and BLs deeper than 700m. Boundary layers 
shallower and deeper than 400m and 700m respectively were not well represented in this 
study. The technique presented in this study was not sufficiently sophisticated to deal 
with complexity inherent in decoupled boundary layers. Therefore, another area for 
further research is the development of a technique to account for decoupling processes 
within the STBL. 

The ability to estimate boundary layer depth from satellite data is a valuable 
additional tool for mapping the coastal battlespace or operations areas. In addition, 
boundary layer depth analyses could provide valuable additional data for local mesoscale 
models that have been deployed recently. This study points to the way for an important 


new qualitative use of satellite data to support Navy operations. 
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APPENDIX A. TABLES 


This appendix serves as a convenient location in order to consolidate the tables 


discussed in the text of this thesis. 
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Table 1: Cloud amount and lapse rate assumptions and how they affect the 
accuracy when compared to the actual height of the STBL. 


Cloud Amount Lapse Rate Linear Y-intercept | Overall Error 
(fraction of BL) CC/Km) Regression | offset (m) Estimate (m) 


75% Cloud 
<400m) -6.5 <400m 
41% Cloud -7.0 >400m 
>400m) 


75% Cloud 
<400m~ -6.5 <400m— 
41% Cloud -7.0 >400m° 
>400m~ 


*_ Denotes "first guess” process to determine the 400m dividing line for the two 
sets of assumptions. 
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Table 2: A comparison of how the accuracy of the technique was affected by 
using Tg for the Tcy vice sounding based Ter. 
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layer) 



















The “satellite equivalent sample” indicates that the same soundings used for the 
satellite-based results were extracted from the full sounding set and run 
separately. 


* - Indicates the assumptions were divided according to the previously known 
(actual) height of the BL. 


** - Denotes first guess” process to determine the 400m dividing line for the 
two sets of assumptions. 
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Table 3: MAST experiment soundings, with their corresponding position, 
time and analysis. 


Sounding Latitude, Anal 
Number Longitude date — UTC 
eu. | rete ie 
GL02 me Dismissed — no cloud 
GL03 aap Dismissed — no cloud 
GL04 ee, Dismissed — no cloud 


GLO0S5 os O8Jun - 0535 | Dismissed ~ no cloud 


122.00°W 


GL06 eee Dismissed — no cloud 
GLO07 on ae Dismissed — no cloud 
ae Dismissed — no cloud 
re 
a 
ark 
me 
aa 
a 
sear 
ca 
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Table 3 continued. 


eat — 
air _— 

— Dismissed — no cloud 
or 
[os | ae 
al 
ar 
ort 
Ee 
3s 
ime 
ae 
mer 
aoe 
mew 
sae 
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Table 3 continued. 


a 
a 
a 
a 
va 
a 
cH 
re 
a 
ve 
a 
= 
a 
aa 
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Table 3 continued. 


aR 
ae 
a 
a 
sre 
“a 
20 
ta 
a 
a 
a 
on 
oa 
ae 
ae 
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Table 3 continued. 


2A 
al 
oar 
ee 
sew 
er 
maw 
es 
mw 
a 
eal 
er 
eae 
aia 
aie 
a 
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Table 4: MAST experiment soundings, corresponding satellite images and 
analysis. 


Corresponding 
palente om Analysis 


Sounding | Launch Time 


Number (date — UTC) 


Sa = nl}; 2327 

GL25 12Jun - 1732 12Jun; n9; 1740 
contamination 

-_ -2020 

_ “2132 
contamination 
contamination 
contamination 
contamination 
contamination 

TelaneeMIAGRBloninen 11; 247) | Vismussedscinus 

contamination 

GL51 17Jun - 1156 17Jun; nll; 1229 
contamination 
contamination 
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Table 4 continued. 


contamination 
GL56 - 1809 n9; 1737 


GL61 - 1150] 21Jun; n11; 1320 


GL64 - 2053 
GL67 22}un 1148 | 22dun; n11; 1308 


GL6s a9; 1711] conducive 


GL78 27Jun — 0241 | 27Jun; n12; 0246 
contamination 


darkened 
ots "2083 


darkened 
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APPENDIX B. FORMULAE AND PROCEDURES 


The following description of the formulae and procedures used in processing the 
sounding data from the MAST experiment was provided by William Syrett of 
Pennsylvania State University. 

The measured parameters were pressure, temperature and relative humidity, with 
winds calculated using the Omega network. The "Omega" winds were averaged over a 
four-minute period by the sounding system software. Raw temperature and humidity data 
were output at 1.5 second intervals, with winds output every 10 seconds. The raw 
temperatures and humidities were cleaned up and then averaged to 5-second intervals, 
while winds Tiss simply interpolated to 5 seconds. 

Derived quantities include height, dewpoint and mixing ratio. The formulas and 


procedures used for mixing ratio, dewpoint and height are (in FORTRAN format): 


Mixing Ratio (w): 
w = 622.0*(e/(p-e)) p = pressure (mb) 
e = vapor pressure (mb) 
where: e = RH*es/100.0 RH= relative humidity 


es= sat. vapor pres. (mb) 


where: es= 6.112*EXP((17.67*T)/(T+243.5)) 
T = temperature (deg C) 


Height (z): 
z(1) = z-1)+(R*TvA/g)*LN(pG-1)/pqa)) R= gas constant (J/kg K) 
g = acc. due to gravity 


where: TvA is the layer-averaged virtual temp; 
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Tv = T*(1.0+0.61*w) ; w in g/g, Tv in degrees Kelvin for height 
computation. 


Dewpoint: 
Td = (243.5*LN(e/6.112))/(17.67-LN(e/6.112)) 


Please note that the relative humidities were adjusted upward, based on 
both observed cloudiness in relation to reported humidities and also on a conversation with 
a Vaisala employee familiar with the humidity sensor on the RS-80 sondes used at the 
sites. The adjustment procedure is similar to that used for ASTEX soundings, except the 
maximum upward adjustment has been increased to 7%, up from 5%. The adjustment 
procedure is described next. 

*** Upward Adjustment of Relative Humidity Values *** 

An initial data quality check was run. A maximum value was determined from 
the accepted relative humidities. If the maximum unadjusted humidity was greater than or 
equal to 65% a correction factor was added to the maximum humidity. This factor ranges 
from 1% if the maximum unadjusted humidity was 65%, to 7% if the highest humidity 
was 93%. Maxima at or above 93% were adjusted to 100%. 

This correction factor (C.F.) was then linearly weighted according to the humidity. 
The equation for humidity adjustment is: 

RH(new) = RH (old) + ((RH(old) - 20.0)/(RH(max) - 20.0)) * C.F. 

There was no adjustment if the humidity was less than 20%. 

Obvious, singular bad data (spikes in humidity and temperature) were removed. 
Unstable lapse rates found in the lowest layers were the result of insufficient exposure of 


the sonde to the ambient environment before launch (primarily early in the campaign) and 
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were thus removed. If the surface temperature had to be lowered more than 0.5 degrees 
Celsius a comment was generated in the summary table. Only the lowest levels were 
adjusted as by 20 seconds into launch the temperature sensor generally had "caught up" 
with the actual (generally dry-adiabatic) profile. 

The general data quality was quite good. Temperature data were eateistentl of 
high quality with the exception of the false super-adiabatic layers. Humidity data were 
generally of good quality, although values were a few percent low when the air was moist. 
The humidity profiles for several of the first approximately 35 soundings appeared 
"unrealistic" in real time. The problem, if it is one, manifested itself as regions of higher 
relative humidity just above the inversion. The resulting mixing ratio profiles appeared 
unrealistic. Specifically, soundings 10, 11, 23, 32, 33 and 34 are suspect- use with 
CAUTION! 

Note: Soundings 10, 11, 23, 32, 33, and 34 were analyzed more closely and 
soundings 32 and 33 were deemed useful for the purposes of this study. 

Low-level winds were unfortunately difficult to obtain. It was often difficult to 
get a "lock", necessary for low-level winds, while the sonde was on board. Winds 


generally began being received when the sonde was above | km in altitude. 


T7 


THIS PAGE INTENTIONALLY LEFT BLANK 


78 


LIST OF REFRENCES 


Brenner, J.R., 1994: Continental Aerosol Effects on Stratocumulus Microphysics During 
MAST 1994, M.S. Thesis, Naval Postgraduate School, Monterey, CA, 60pp. 


Eloranta, E.E., J.M. King and J.A. Weinman, 1975: The determination of wind speed in 
the boundary layer by monostatic lidar. Journal of Applied Meteorology, 14, 1485 
- 1489. 


Davidson, K.L., 2000: Personal conversation between K.L. Davidson Professor Naval 
Postgraduate School, Monterey, CA, and the author, 15 May 2000 


Durkee, P. A., K. J. Noone and R. T. Bluth, 2000: The Monterey Area Ship Track 
(MAST) Experiment. J. Atmos. Sci., In press. 


Gaynor, J.E., and P.A. Mandics, 1978: Analysis of the tropical marine boundary layer 
during GATE using acoustic sounder data. Monthly Weather Review, 106, 223 — 
282: 


Hasler, A.F., Strong, J.. Woodward, R.H., and Pierce, H. (1991), Automatic analysis of 
stereoscopic satellite image pairs for determination of cloud-top height and 
structure. J. Appl. Meteorol. 30, 257-281. 


Hsieh, W., “Atmospheric Convection”, 
[http://www.science.ubc.ca/~ocgy308/chap9/ch9_c.htmlj, Aug 1997 


Kidder, S.Q., and Vonder Haar, T.H., 1995: Satellite Meteorology an Introduction. 
Academic Press, San Diego. 466pp. 


Klein, S.A., and Hartmann, D.L., 1993: The seasonal cycle of low stratiform clouds. J. 
Climate, 8, 1587-1606. 


Kren, J.R., 1987: Estimation of Marine Boundary Layer Depth and Relative Humidity 
with Multispectral Satellite Measurements. M.S. Thesis, Naval Postgraduate 
School, Monterey, CA, 73pp. 


Kuze, A., and Chance, K.V., 1994: Analysis of cloud top height and cloud coverage from 
satellites using the O,A and B bands. J/. Geophys. Res. 99:14, 481-14, 492. 


Ortenburger, L.N., Lawson, S.B., Patterson, B.J., 1985: Radiosonde Data Analysis I, 
Western Division GTE Government Systems Corporation, 1985. 


79 


Rogers, L.T., 1999: “Refractivity from clutter (RFC),” Presentation presented at the Naval 
Research Laboratory, Monterey, CA, 4 February 1999. 


Rosenthal, J., Helvey, R., Battalino, T., Eddinton, L., Fisk, C., and Greiman, P., Predicting 
the EM/EO Environment from Satellite, Synoptic and In-Situ Data Sources, 
Proceedings: Electromagnetic/Electro-Optics Prediction Requirements & 
Products Symposium. pp73-86, 1997 


Rosenthal, J.S., 2000: Personal conversation between J.S. Rosenthal Naval Warfare 
Center, Point Mugu, CA, and the author, 13 June 2000 


Shenk, W.E., Holub, R.J., and Neff, R.A. (1975), Stereographic cloud analysis from a 
Apollo photographs over a cold front. Bull. Am. Meteorol. Soc. 56, 4-16. 


Simpson, J.J., McIntire, T., Zhonghai, J., Stitt, J.R., 2000: Improved cloud top height 
retrieval under arbitrary viewing and illumination conditions using AVHRR data. 
Remote Sensing of Environment, 72, 95-110. 


Stull, R.B., 1988: An Introduction to Boundary Layer Meteorology. Kluwer Academic 
Publishers, Dordrecht. 666pp. 


Tjernstrom, M., “Marine Stratocumulus” [http://www.met.uu.se/eng/forsk/marine. html], 
Feb 1998 


Trehubenko, E.J., 1994: Shiptracks in the Californian Stratus Region: Dependency on 


Marine Atmospheric Boundary Layer Depth. M.S. Thesis, Naval Postgraduate 
School, Monterey, CA, 88pp. 


80 


INITIAL DISTRIBUTION LIST 


Defense Technical Information Center .....................ccccccccccccecccccceccccccecssccesceceeceeeeces 
8725 John J. Kingman Road, Suite 0944 
Ft. Belvoir, VA 22060-6218 


BSB Kem Ox Paint tiny career eee ase cc cctay cues scotsssousessecsessseveenseswessedeesecseceedss 
Naval Postgraduate School 

411 Dyer Road 

Monterey, CA 93943-5101 


INVA eaawalinn geo lite eT) a ees eo on oes ov aces acetelnccinseccscecectecdeselescdeesstecedeavedesteee¥e 
1200 Greenbriar Ln. 
Arlington, TX 76013 


RA CONT Sey pa e711): JIE eases eles or ae So 
Department of Meteorology (Code MR/De) 

Naval Postgraduate School 

Monterey, CA 93943-5101 


PUOTCSS OI alten emis ACM resin yee soo reece, ss <sau5 ebbeeeshee suns sassvonsaceeeeebadecones 
Department of Meteorology (Code MR/Wx) : 
Naval Postgraduate School 

Monterey, CA 93943-5101 


DAA Tey 7 Te cco ese ae ee Rc oe la 0955 50s SOW ahaa gesnce cov cancessessseasdeeiuseeeeseets 
Department of Meteorology (Code MR/Jr) 

Naval Postgraduate School 

Monterey, CA 93943-5101 


EASA AN OSC IMI calls eee NN tee nets, eS dnd ade aucun ncaa sabe vod evensneacgenneres eaeet 
Geophysics Branch (Code 521400E) 

Naval Air Warfare Center 

Point Mugu, CA 93042 


ENG) le reg Pe VON 55 a5 ore NER elas bog ng’ nt sla <a odadhat. ache canchideies iceeeenyesegsee secs 
Geophysics Branch (Code 521400E) 

Naval Air Warfare Center 

Point Mugu, CA 93042 


81 


= 
[> —_ > > ——- 

ee = 
oe = — 

—_ 

CT — a 
= 
& 
= 
_ a S> 








60 vm Cedb ¢ 


6/02 99607. 
















































































































° - 
ee. e “e — = « 
. - e. - ° 
- Sal eee - f= a: - « . - 
- et w - the ee ee - - e 
a ° o es w~e e “ewe - => a = - ee 
Py . a. al - -° - -_<- ° =a 
7 ° ep - *. e ° one ° 7c « - - 
e e e ye ° ° ee e = . - « - ce a 
e = - bed bl - - e - ee ° - 
e . < 2 * es oie r] « - Py = a ° ans 
= s . . . = . e e - e - e - - a, - e 
. &e . “a we ° e e e .* ~ « « - ce - - 
. ‘a ae e e e = = -- e Pe ~ = 
al e a . e - = *£e a . e “en - - ee 
oe, * . . °. . . a * oo “@ 6 a « ee re = <= 8) se; Ben. -ae 
> DB . 4 ° ° « ae » a e. - - °. = © « ° 
- =e . . pene we i a) = v ee . . ° - e ° ° . - a -- = =e . - ° 
“~ . - = - a . os ° . . . e -” ‘os ° . - . ° van ey eae —s.% « - * e 
- 2 . = - mF ° 8 7 . - - 8 s . ° . 2 .e ° e . ° fe. = - en a er a) . s - 
o Be ~ = s 7s Cs bed e- . * . a il - ? ee n~ e = ee - -« . a 6S 6S ° at eae - - = £e - s -« ° ° - 
= me Ley a ° - e- ae os “eae aa 8 « . « es . - ° < <5) =. a e - ‘Vania - e as = ae 
_.. "@ ave « ae - = bi 9 - ae 7 = - = koa I . e . . = _- e « a . e e — “ee. - “aww, te - img it 
- . = = ~~ e == . la eae . bal - bd . 2 Li . had “ - . e « on6 es e eo mie = =» e eee - ° ° x =P am Sin a e ate - « 
=e =e -« - . - - ae “ « ° . - a = - a - *~.. - . . = . e « ee. we Ps . He - ee @ . ° e — eee a . ea . e “7 « 
== * w= == “=. - Sd ~ os = = ee! ee . - e . oF e Ad = . =—~ - ce a ’ de - ° ° 
« -- wae - - - - = e-« z= - ™ bad - Ya - “. e - « . . - we ae f ewes -@ w- -< t= 
= - = ec. =o ey ° 7 2 ° ad ° ave Py . e - ee ° e ° - - o ame 
- ° - om - ‘a “=<« se es ° * - ~ Mi ae o% bod « * . - ° e - » . i e Ca - a e “= . - wom e . - - Pa ° ~ eo 
* = Sede Fee rerry Lbhalle Sune Ema ee ~ - . = aes be - : ree es ° = ° * - - . Z - “ys. ° ° . ° = 2 .° - a © “ 
. - -° Care - « te ee * - + 2; = -- = ° be ae - of ° - - mes peice e . a « "e ° ° « = ° @ee - = e - © * om 
= « - ee. - =. = ore - - - 7 - = * ° e - ° e ° . ° o - °. ow »4 o e >=. = fsa’ — Tee Sane a - -@ 
*- o- om . “= = rs’ = ~ . ° - +e . . e. ier %. - ° -« = - ° ee “=< ¢ ee Om oe. er -« 
- - - - e -- B25 ~ i: = we - bold - = . ee ° « . « oa ae e ¢« ° a -=s e “ye -“ ee 5 on « e- eon a 
- - ~~ = -« s “se «= 7 = . . . « a - ae = e « a ee ~~ = "=" = co ome 
is r, - Oe r ~a. . ° ee = . ae = a ener ras ma wee eS < = 2 re eae on eg 
ws -° ° - ~ =. . . ° @ e « - ‘ = - ewe a «- ° . =? we ae e On ee een e - . ~~ oo 
x eum a > $ e = ~ * . - 7 8) * on eo - - - 2 =. = os “© we ad - oe om 
ws . a - - ” - . ce . . Zz e . ~ "« - ee, ° « = - « wae . om « - ee @ s -= oc © ame 
ave 2 7 . - ee MY at ° . * . . - . eo - € -™=- « - - Ft a Swe aa "0 © amy e x. 
A 2 . = - + ° - « = . - ~ ° 4 .- . - ee - - oO etree a. ey ~ 2c mee o waa - e« 
e- en aed = > m 2 is, ba se Se a aie . *. © ar) 7 - - -- “. ma es . =. <2@ e +2 e@ @n- 7 ane 
é ; = ° - ="ee e « - . ° a . & a aera me ace << ae s o.0e ef « = e - — © ecm ° ~ cm e 
oF tary - ° . . ° “ © « ° e Pes ane* -~ e@ lee Pe ~ 2 es ae RM ie cs ° ~ 
= we . . - - Fi ee e- ~ ee . ee - « . os = a -* a oy . et ee ‘a =ee es - <0 = eee ae oa ow "e Fo Ow 
eee er wet = ae = "+e « = os) * e- . - e eee a ea - - en ee « ae 6 » PONS line. a aoe a er — en 
- . . be * =e 7 sae. - - - _ . . . o~ - PnP en ®e e . ° ° ° Ae Sr ee an wenn ss ae. AD om wy 
8, Weee a ee * “ee =e ms . = a « =" .@ ons e oo * - Ear o wt | eee -«. ° "2 e 4» © wea o- — eet on ee ° = 
e - we ae - - - =.«- Co] ~ . . bas . e- « ee - = e Py . e “es 6 us ee ° Ps °. Peee ov.- Cem ow 2 * .e-m, 
“= wos gn e--. ._ ” di en tee co] - os Ly - - . ° 3 ww e « ~ - =. SS eee we a =#5 el ee cL ae sy fe « 
,. 06 a ae = ’ . oo Pry - . ee a ae « . . “~ ane, ay aes e . me . as © ea Sel Tie ae Sea) eee Mea ome 
< - - o. =". . = . e® = oe Be tess = ae » . id ’ o ee “eo 2 8 te 0 8 6.68 wan . “eo 80 -e e ee . . Wee Oe a) 1S wea an a 7 On ew 
- = = -<- - - ° = . . ° a af @ ee ° - & e = ° oe * 22 06 ° “=. . ow - se 94 ee w - e "ee oF How ee A oe DP oe ae 
ia - a - =e. - 2:5) =e ec + . "= 0 - - - - . . ° e “J ee Se) .o on a . eo & “2 oy eee . - *eeea.e oO, ote —_~ = meen pe Pt 
o a ~ = ==. . eee we wre. & ed = bd baw! tars! v= 7 e Lear “ . * ee es - = Son, -26 =e . = we 7 Wes oe ~~ meee °c wars "Bet Bee e 
A mee se . “. 7 “= 8 8, es: ae * e " e- on « = « . ten e ' - - ne e. Cama co 2 - Fok we ° D2© oem Meee ne - @ 85 aw = wa =e 
- ° . ? ~ . ~ « = * ee i . ~~ « ° . - e. a es Me 8 8 er on, = 5 ° e ee eo - 2O@e “ee . “— = e's @ 06 « sco rae 
Pao a e . - ote a = =. o e -ee . e -* ee -~ s e . - ws ¢ a = oe 2" oy *. ” =" o= . -- te =. Fteme ° £2 ©, ~ wee e - we 
"2 « =« fu @on cose fae - - -. - rn * a . . «@ ee ene, e Se-n 2 we - 2 wee ee Sito = faces ~ Se ee om. - aie.) 
etiam: "= awe - ™ . ditt: = bd - ° . ~ . . ° se ot e sPeezae wo , ue i ‘ae ‘oe er Ow, eo e - 2, wes 
CC Ye “= 6 ee ¢ we id , bt od) ° td "- - . * . . ate One “e e . ee oe ‘ae 5 Crs 6) gale 8* epee ore WP oe Co ome am, we. Se Pee ORAL on 
2a “= 6 ai a St - ' - aie Lo -¢ ig . Soe Me, ote a 68) os Se a) eee a a e- = me S © enews .. a Come 2 © - ¢@ . 22s 
—— os - =. . . = 8 " 8 we bed . b ° a - ve e » -* aoa or Bel . a ms on ee of cc pare oo mt le Se Oe Oe en "2 fw mona ) 
aes a=, - =. 7 - . . ee oe “mae . - = © Sle « e ane See ae e Pnae ~ “Se = © 1s _- aA ~~, eee @uPe - “OTe 2 wae 
=o * 5 . . bd . . =enna ce es be « Sie hal e - .* . . Ciao ow . we « - ee: 16 e e « es . - ° - wee @7eee © S.2me ean ad tse eur naw 8 oom, 
“ee ae . so ue = oo on ole! 7" e4 is mer eal i 2 * - = “4 ae 7 oe ’ =. e Vine - e#72e Pee Oe? ~~ 2 = ow ae ~e- ate) e oe -— -2 « 
== = 2 . - se = ’ e bs cf Sel ve LP . - ” . WS cee as ° - e * -« =e * 8.0 6 we °° w« =) 9 @*yhrome = 7H 1am e wns ee e _ 
Pay mime = =e om Sy ite ee Sere) s as Lat ea bed - .- e - _-™ . ~e ot @ - oo w ae © PSewe wee . me, i ue « wee. oie) ee “Om we 
is = @ - . Sy Ow =. « = ae *. e . « -? . a e-. os . . e - a =o ene Ommy wee © 2 & *, ~tgr-. Oe evn, *. xs we Xm © —- am, 
ee ial ws) YE Me é on . = a) - LT ae = ou ea ° ° ar ée - =) Wea ye we eee oc =. . 2 os one ~ Fi a Fe wey 
= ewece “2. = oe be Ue -- ' Caked ° = « . =e . Sie: > one 8 -. ee ny = & e =. “oe ou ee & & Ore © we em. 22 MA Amys as tome me 
aM a . El e = “¢ - . . seri a San x rd) es "- . . « eo « . e ~~ ee © © ane - we mm, & « ee NwTest ame 6 biel le on cerns, - 
ete) Seced oe 7 se C1 ere m lead ae - 8 2) od ee © 08) ns a eens - ear © e & %e a LY - se 2 Cehs seme co w © ° me ow ae FP Ae: & OPAC oe ws Se, = gl vane 
ome = ™~ . -. -—- °. ae - pul amet | - . 248 -< - ° e ma au e - ih er = @e fe ea oi. —- ~*% @ Soe) ORS ne ean o Crm “sas 
Peed! mel eas -« i) Ltt “tne - «© a | . 8 ad . a =e @ . ° « i FP es- oe PRC 6 my dae 28 oe meau-e ore ae “ee oon, . ~—™® ee wm 
“15. ¢ = 0: ome . “Sree eo. s . pitta LU TE # - 6 Oi a esa! . . . “e ~ 2 2° em oe . * 2 oe 6 -” ny « 7. Sa5? | ial Om Pe sam o == & 
. = ~ "Pw ane - * = re eet Se uee o te ms, See Sieg J e S te « . * ° awe a: oe , * e a er e - -. a we a et 7 ee an oe ZHene 
=e. * "ee = 25 = oo = . * ° ° cS) ieee) on ° . . 2° 8 % a ee . conn mh = noe °8 wary ec bd “@ © ee ee x ey Gs om we iw “2 ym -~2 pwr 
a Foun es as = "= "2s =~. © +s aor Se eee ie ae i ee _Y . . . "=. ° - et2e mee ° COM uF sey . . Sa ce? wits — awe « id Oe Gt om ww © wes ta 2 owes 
er a awa =e. wae -"=em e w= ° ~-~2*e ° OS aa toy, - } a Ts, aN a leltays eine se a wan eo". » ® 2 00 or ony "Fem. e~ 2 w- wa a Pm. mee y 
. on wen omnes fe "5 =. ° . ‘se Li ~ eR ag *? <a ay am - e ee et ee 2? ae = . se - om ~8. 6, °-2 & =e a or PRM egn & BAe tean -~ 
ce =. eo =: sees - oe ow b O . - = @ et nae o 8 #u “se eur e ° - ‘ee . = aon = - 8. Ponte | ee . 2.- “80 eo Om OOPS Be 400. Pa Bad 
tee. a . e- eo om te “. “« ke td - . or “er. = «, es ° ve se oe ob So @.y- Ow <a e eet am eae in > - dt te CO fe we ag 
a SP0 As omateace = - e"= we ee ° ro ane = at SIP LL ~ awe a . . Cd & e ~ ee few a wm ye %& Me a =. =e estas tos 22.0 = PRO 90 Mee we ot °= 2@ 
fem y eee fg ek S = et 8 o fen ° bate is a ie atlas . ° . 4 o- “2 ope o Pe ews a ee OOO Oe fan oo ec © eam Beer or Cae e° + wees wr se 
a Sia (ay. ue pte « - = ea - es -- . aS . . ec fe. ‘ e e oe om fe-e oa a a a6 os ove te 8 8 =~ e@ee Fe ete PG cal TE tg - — . 
= an «ta ~ ". we wee -.. ot di ore) Sy aoe So) (8a = bcd 42° « 7° - » e@ . ee ° « wa « e—=e a em. etn = SF 6 me Ld 2 = © age - om = a Lee gee 
nme 9g. ot “Se amen ae 8 eer fee wry ‘ " ee ae ~ > md *. * - = «6 ° . ° > none e-@ Sate ghee 6 ue ae om FP eo Soke wm SON SS See ie ae ack - « - ween Lf — me ee 
. = =. = & = - sree ® - . *e = = Seem Te a - © ee 7h Ong « . Ce se . «= = . e o2a “= @.=e@ On oat « we -=$%ee, eso we te a 2 eee eke ee 7 es =m oe Swe ome 
ae 7 me eg a a “ade ene og ° a e ae Te eee: of to oar - é - ©. e hers = con mate ete ee Co . . © wate entra ~e “°° Oo meee A © 
eae mina Ca e “ast mee picid a ORR Sm al wes ele ee ad * . . © O. mHe = 8 ow ee 1 SS Cia o0s es e » or See OF ab 0 ogee oF Pete eo -— we -- fa Siet-@r. ame a, Bee ad 
oe oe ae 7 27 SONS AP Cpe ey . ee. - = as = . « en. a oon Capea « . Py © & wae rm ° eer. , Sot ee, fee © om ea "so a Bee © - © cope ae 
Resi env oph ian = Pe. ae - © = 2 wane = Sm 85 a ne bh I) + Bee (eee A = . ed = mo te te fee “Been © Gye en Ge ielalan me 8a, 68 28, mom. of “= we oP Be ccs -e ear 
Poo ot =“. e. Lit bet wat Seer res oo er @9 ay igel =e - ed Oe areal eign ap ee ay ale “. aa - = one ar ai etenie vee e % *.ve= © ee © © Omen wee of & Se) eee oe met, 2 ome i= 
ee So SO MORN yaiah ates a F "2 OF msee® Q= on" are - Ore cntied *- aaa, he o. ee we . © Oneser “Ra@ge a eee w Bue, * 8 6s capes me Ome = Se hen ~ wete, = 2, ee Wet e «wes 
were, eae fee es wae = ° Be © «me “se ~ bi Mee St se wee oe - . *. = . . 7 - = tee we, eet - =e ce Me mee Ss Comune cue a) s% 7 * ee on “es tte 
Ca ee ete, "Pe. aan a a) 08) 6 *. me ows . ° or Fee “e - hey, Py = ee ‘ea . = ‘eres a a Boere « . fa ane? mew « - «= a, He? A, — “Pew © one ~e “ma . arg. “wae Cus 
an "etm me & BL - = 7 "ere te oa os . a og oe ° » Ps = <i Sial em oe a =“2e.s, « Cee . oe eS eee Tm em 8 we te fetes wr em em ee, oh ee 
Say Wages = 6 -. - < Nt 00a, OLe) see ae ee ° oe 2 a Sie . «= « - e wngae 6, | WP eSe ele e a, en a ewe oes om 8 Aome oa wee os » A Ore, regis = oe - 
“6. “=. a ~ Be tree 2 ee = 7 2-0 we “ee cae = o- ee -~¢ eonn e BUA SP gs cis ae few sme ee i oe Ss ade wo" Ott fe= wee owse 2AM An © « 8 ay Se Pee ten cee a @.« =~ fe 
- ide tt tier ry - . ede .- 2 «68 = a e- s O88 00 oe wey WR isieer ae -_ we See *, aes “" was > ——e & » OPa® 8 lems eo - SS an eee 2c "2 es, - “A ow ee 
CO q Pare Sus “a, eo. “Oeneeee 9 ounm § iniealet Jane tee as Sata ) - om 8 SuAe=n © ag, *. =e e . e sere er . mt he ° “290 . Miepegae » = = bene Ln) So meee Pam OD pees © heaiel an] =* we a Son ee ae = 
— BEN CEE tee alacant » seeig Shae ree: PO astalay ee = 94 ars! . eae - L ee f . * 6 rie ae e "ese 6. * = ot . = © ots omy, © #S © Gen, ate Sk ee = ° “Stee 6 = ac 
———EEE . pire Maney atte t gay ene. = = ete. * «0 ° . me . ° %. a" Fe acy PxT) ” we « _— Us SH we eaten aed pot Pe 0s Mee oy ee Rene ~~ a --e Cee Ac® malo cgmen in ete 
———_= co ee «Le = «= = SMO a a. ies eae ™ a 81S a Qe "ow oe . eo ee oO cents one yg ° oe Ley “g*etee “= oe Swe w ot space = 8 et trae tp he o-905 ~™* a we we 
 ———== SaLaee RSPR R@all oma g Tes We =O wee og, + CM ar s . * i a ee On we era a) spe re ee Maser 3 A .3 . 9 te eS seh ee meets § on OO We ee om om —_ See Oe Shei co °We- «6 AORO wee ot 
= —== age aaa Fe taeu ren "e ee . Len) Bere “Sh Oren « ° ww Fe © ene oe ‘ aa eae, ies o cume . aie “e: wire .en fee oo 8 ome, ae nce «OS eee Aimee 8 Aa"O% Roe, © yp PD wee 
oc on) aur oen-e . ie Sw. © euP ee La} ce om ae bth dal = r= a = “<8 =4 © . 8 eee a) oni y = ce RaW e ws 8! dae ee e = ~ =" Sel ea, ee be anttkel alte PY WM ee ws *@ se .@ 
<= — fe enamnienese “ee _ - - Pe sen ne net ue boat ea =" 20 %& - es ee me ® na : ome one erie =e a © wen, * watee “sue et te we Lage Staten « 
=— ™ mete Me yee one, S00 F ene me met cone a S68) wa 0, ole e we “Om 0 me oh See See oem, Py sale 4 Pra ' = . . vo 2 be Bow wes cmp Pe emeen, B FeO enortn wa oo om OO. age ye on ao, mace 
Cc —_ ——— ES Rana 0 ela) fepanasly hchder tlt tele We - ° on i . te er ir tee A te tay a ae Pe eee = e SE iGieee in. ce “sem wee =e wet in we OPA ws 
ao CJ os. ap ciwes ee mien ee at es Sie! 658 ening S = wes ware a oe » -e - . . e z 7 ce ‘ We! “ex sigavecsete we Penecsey aucmiallge °. 6 
co —_——————— wee. ath od Le er eae I Ley ° © ewes see & . = BAG Fe ele) iy ". See NC a ee 5 @%ea~ . o+> 8 * Ree ee as ee er =e ( 
—_— oa hele. me = 2900 Bee oy PCs ge 5 Con meebns nce oe 6: wwe 2 ten be "ein a ° ° -" 8 - 000 oc am « . ow *, * tee 8 atue © o & = 
a5 © : “2 "8 -y ae Se: 6 Aa0 ee ac e=0 ee estan &  ebne “== “2 of <0 esevg ue * daw . . . = "e 4 i . - = . a as Gane ° & © see 
——.——— p aerate: sagen (Bae gee. ee dee Paha WS Wt Oe ser & fe an ty x16 Pe el ea cde eel US (fa So Caer Se a) - 6 ee “e ° Toe, “gee Pee ee oe bd 
4 — bs "@ © wm ae PMs 0. ona, Seis) 089 seam ce et es av 6 oo * sw 8 8 = 08 oe ae 45 - ae .e sere » e ° - © = enn, ae e : . Oe tao htt ee 
= Wy, Or ge sseslaer@an, as Let baat sce we «08 a5 0 te © meses "ao we mek one Fe fone . - : ote ° see ne. * eo fF ae, ° Ph ° “a5 8 ’ ‘ Ont ar eer - ee . ‘ 
O a © “s Se ee 978 =O Seyec cm Ben Pi ae aay “2 Gree abs s 4 We teens oy on sMer te va oon Se Colette ee & =005% 6 eree * ‘ . e ¢ = ef ~ Cy a! La. cat ef a 
ee ~ “se e Lar "ON eee ce wa” S, = a 6, 0, a vat he ° : 2 6 65 Bas OL ary tes a9: ise bee rr Be * as ry « ie or ow Caray aie! ae . ee ene Leper o "Ot. ew ot @ 
rr ee s)Ge Ute Oe!) wel game _ SPP t= es oh ge: men = Oa tare og eM A et ww, ° e . n~ Lota 2 jee "e e cose yg . . wo oF © Bees oy . = ces .. rs ene) cel te Wwe proce we Pa? =e ne, © Reet © 2 werrie —* o-aaied- ate 
F z =— — © =eten~ le be er Seale, “Fert ow os tem « Bees Gon @ "sf 2 an 8 - BF gas co he = 6 = % e a wie an d “a e ty . -- e- 5 he Bs as ° sire ot. of sap ven a a Preeti ae e ieceen “wheteires Ome rn. Lord ek =. Gon 
CS Cire, ee eae rag he Sh ROAD Oe Bite on “4 Poomrmy 1 om ee a ante oe @ Ge ". " en 6 | cs. . aie. ar 2 rs Ps ICE ee MP8 eee Ne age tl ee rn no, ane ~~ eee 
==8% a2 "=, om tye ote © na =/ ese © 1s ot "Otome wut wae © Ne wésaaue @e ee oe e ~ @ - . leas = 8 ie * ey @ e ns fas a cane Se, ee % 6 mM wee SB 8 mage a « @@ ~o% Se WE ow wang - Se oa. = Soret etlase Pern neon ce +24 ~ ome ™“~* 4emes 
~~ = [0 @) my =aeh © Sen © we pate pale o One ae (én . a % ‘a eae e ewe oo. e “~ _ ae e * OMe o S ew os 2 « #e* oe © oe 6 us fo © .@ Fee, ing a Oo ty boh-ebe inate nt © £088, owas ety Foden Sige te SO" Baus ob. oe Sad “tune one 
——_ . AG% Srey Se, “Esa ote ae at ate M *: = "ee 8 6 ee r ee =" @©8e © ww ©. he tw e ~-6¢ ROS (0t6's iw a2 4 « oy « e oe Pameesa 8 ~~ © & wwe ag, “@ =OMaee 208 a, es non Sent. © Ccteqee we, ? wes Motes, OO gay ram’ pt ome rem een @ = wet 
YW a hen Oraye, seas ene @ 4 ie le an | Pw Peo "Sw . ww © «4 ©@ Pa - ee on eo? en een & ya os = Sie ae. ee «© wet eo omen . ww 8 «me Oh ne “tangs y (OMe me we, eet Ae eume VO = eee rE we Age Ort OR teen ocean 
—_ eae ae Waseem e woe Cpt SOrtw eanrore gy bcd DB Soaked Y Sa ® . os - * he ey | ewe @ ° ~ ae Brewe we gon ee ee ose . re - ¢ =. @ %85, a on. -neriie *, 5 fee Pare 0S Mae wees eo ce ae PP nale iene AO eee @ Com, 
aye Seu = wom. ° . s - ° ° - > oer 6 a. ate « * 8 ee 9s we Mne nee & “b= wie eke se wee a* ae Sonetesce @ e. 6s Ses 2. Ceo nw an Mante® ©e cae Se & mee OMENS 6 ocmem tw sche ae ee WM. wee 
OQ {™ = Utne om, ae . iP Wig. ae o8- oy “ee 8 e = ee ate -. = s Aner a ives = oe S20 8 ey “e Pee dene ts a8 80 ase, fen = Sas chee na* ter « = Sane ae) ot * Dera, sen ga ete ee ee. 454 
Bote ceetin pa St aen =. “es * Se- «6 a “ @ . « lt apt Mees Se a) ha Sarg eee © «ec fe ae ten Sf 8 epeg 8 tree yg Cot whe gece Cd ees OMe ae ee age Sere at et rerew ss eosenewen ee te aan a ean em 
= NJ ower e "Se Casings s = e6e.8.8. = ree D a ee . arene ° 26 seis: ®n0 @ne « > 4 - e ee mew ~ ° ote eve e eee “° Peon $8ehe at og Be om - © © 8 ow ~O-w tele te See of o, = POS oe ye NOB e net On fo + «0 ate? MG as <S° Po mem © as TPS Cee © 
Q EEE gee FOCe oon mud ee eee ‘eas! a ce @ » . ue © wae . ee me gw at = oo e 2 - & . o a  svotee oe eda Le ee a @ oF oe% con once om © tee % eV We OV POO Mee a thee hee erent epee 
= te yee ge a are eee Wes. ie wit athe o. ” foal iecim tg ae : ea eee ewes © es 280! S wenn Ao Ces boasan 28 Ht. eee. LP ~ A? Progen Serevent Se Kiel ee tee en me 
- a — a i ee eet et ee ae ow ow ° “Om 6 Oe ents, « wee oy, = ere, 5 Ee ee Os it a pee <8 Qaere gt, o oMee- mae Swe wey oeny =e ne ns tie Be 
ye —_-— co *Were eae sen > we ete a er Ct ones o a eae) eh ate ne meee MieWe ce Seco me tes 8s o8 . Se .8 = = © 2 ce, *s e ee ce oe eae Neve e ree ett 28 gon, Rotetend inde 
ae) — oe ee ow = S2¢ © ® ’ o me e co ew * Sie ere ea me. Rot Besta goconied BicTeR we. ne oN Pores seen “maneces Soak wee -~mee SSS 
ies, ——__——— ¥ aM wee oe ° e-«6 OF) chs ip ate: eo ane 2 pee irs 7 a eo2e it Penis” ec ote be se w ae we a mabye en ee 
»- , — ice i . ee. - a oa ite bs © we ® 4 Senne Py ater, ee . Ponce Lie ee ‘i sor n — o ets tve wa at. An, ou 4 © we. ° “tap PRs wo pan Pee en Om Ot Oncmepe &o.w 
-——$—<— 4 -« > mw ee te “en. 6 “@ © We « ~ © ve vu fastian a oe Seu Py eet be oe Bees eae te wen oy ; . wt eer e ve -e, eo” Fae Oot eo eg WOee De aunne 
= = __ ee . 26 Serum ewe wt we Pe sf ek tpn Nre og oe Rae “ee ie *eeeq eae . 4 : a | “ aug p = ae - ~ wh By Fe ears pe San en FO me oes oe oh 
—— & =a Se OS sca eee - aaiteg Sietiehete Meese a. Paes eg bps Wesiacrar gee a's ay iemiaes ses . “eee cee nate shane * ie Rete we on. bey att Pts 
— “OA use Paw OC ne Og #%,0 ee pene o ° oe 7 “" e.0 tee w e Beawnye “~ @ we 8 « ° toe pep Spay a Gow o ees <e Sito pone, ert RP mer ott Tommienen 
SSP Ps ©. eee e * - % eF ee oe. t Croum « Se 6 ee «ste MP OreP. oe. ae Set eee ts 
° ete opens SLTRe) 8\ ie earae * ta 5 thse Souler auecepe cc nap hy : ae Pa Palette o * ieey tev De ACO een ett Are eee ot. se e-te-aptrdenan soe SOP S00 stem Aare agen, esl me 
enamel iat bac) see ensae, (66 en gs - Salen Ee aie ae Cima GOW Ae OA as Se sa etoee oe wens eo Oe e . Netw Le snae a Aten ses PM ate oot ON wane, SS. SP WO w, * Sena e ee lewwan gone 
pe OL ea mbar t-| - to Oe) aia ie -° wee wt oo 6 tom, oF es Le eesgis 8 stam sm aS i: : 8 Or Sema wes 
" @me eo. se =@8 on. ese e = 6 & «es Crary ww ee See He meme one Pe 2 Pace wee Meewe" a, Sane 4 oe a Oe Th Sin ahS) ree ern e*enee af ’ : - 
WSs fee =. *" Reman e 3 = We one a ge *, *. “ eee . oo "so wm ° @ 68808 “ee eran, en © os wat 4 ee 5 em 000+ ay fe merey stun o en Sane a “PEO em erst 0 gm n pene os 
drtth ed To or ge " See OC Sir Pee Wid ers EW ese . . ’ Sere = me Poe eRe oneness any “be Seem ge %&s ~e 28 we Pe ors aed 
of ue wy be fiit.e iss, Dace eee SOeW Mg oy iY stacee * = o ee OO oe mae pM kteTy ~ Oe we me a ie ow cde SOc ae as 2 op my ‘ Be ime, cote uge oie # eae, 
FA ROEM Wrenn e bad ry BOE Cm Gon, den rete one oe ee ces See DOIN iSea abe oe - oe as ° 2 - ‘ene “C.a8 eo 8 © en ge “mates seas ar) * fa 
£05 whee © Omen BEAK, gs Aor. Oe Che rr ae om aNot elect ee * ove a ert asians Pres Senlece e 
ue oe co = = we. Me @e. . — . e ae ~ ~ S, 
= © © erase ei eet le neha ut ney wem oe, Some, ee ese ow os ete SL et. a 
=e6, ate Ore tev olen os © ee tne, Mare we te - e > os ** % eaye e , - 5 e 
-" eae heme RL oe . SM ve We es © oo oe. ae . e *o Ls ontue ns kee OY, “ 
“we oe ey oes Oeney tent CC = om wderemte . Pore Ln © 20 oe 88 6 be Se eee , © e te - . 
Siew eae a tae, hia Hag ee te EE Ow S68 O’ aive oc - e= eee oo © fae 
ee. 7 thee oon, Ct rs “¢ ° > . 0% ° Se Sr 
o ob aye har i Te wey ae. a aniite Sreee ote oe oo - ® ~)8 WIE, Os ag 
Fea eee seth in x 790r Flin, NO me - Se" © @& eke me . “= %eeesn % anh, . 
Nf -e.e as © aa. omg Nore LP) oe eteem tam , oT N08 OF men “P%e  te oe Sen os ~ye “en = 2 Bua 
“% @: 5 on, eo "RM mro® % = ae se tm Pe ems eg ea ea Cr 50) en ~ &w -‘s 5 
fas OB op pom *ee “~ ewe Peo ene « ° @e e — * sen 277% wesses bd oe.8 of . 
ranee See ee ort 1 ao? See e “2 evens 3 OR meng e Ome Oe- 8 ete © Mc i 
’ 02 ww te te “ey «¢ ome ong, & 8) ie ce totem © ore w er) pe Wes we Pang, = (las ewe om 
mee . of umber o . e . «6 eet Oe . ohm 6 ate Weises "ete psy « TW ne ee 
Fe net oeee : . rei Tey eee . -e er eaten, Fie ene LeO a ewag cas sce SeSAe 2 O88 eo oe, ae a 
: oF eg Oger, P=8e0. 6. of at, - A we oe eo 2 ew soe neha % We Cal as PS ieee "82 lew 8 aime 
Sr eet an = 2h Og . o aSury ou , Oh ry es + w= are * enee., terse Se © b& aeoN en an ow tee er Ty) ie oe a Bang 
een Gu0? My p0nes as, , ? Pare wee aes a « Pee e eee @ ene, OBehe + wage ofPe « F c@ 
o* me's th SIP 8% eS OO Moon an, Owe Oe ee sms tical A RY ol i ary ee ~: led STS ee Pt Te a) s 2 OM. he 63 oy ve hon 
8a Fant eats Ned = at es se ee we ace ts 7) ws te es =o Genome UW Sie" ‘we 8 ths 
VOlEe a see mye vuee Nt?! Neate, ett Be * Fem nw i ee "ae %  *e 
“he es " ¢s © ope Pr seme “ee ae e a8 =e & =v. e ” e & wor, es ef 
ee eee os on Ob ay ene + Yoree Fare ° em a on 8 et ewe oe . Mr 48 6 an ea © 36 pwn tepte ° way? ww, 
Bom =89Bee wy O Pee ute, te MOO PA0 Ol ccccae coe we ef ~ wet Bh eae enhaps ej e tine. « on hdl eet a ah 
het eh  L Bewo® sy . sre t8tBPee Le » 208 WON Neuse Set eile Per me fen oe eee et ., uw oc O,8 @ 28 on Pl, a ape 
es? 6 Pain e poae one COer ome bm ole, oo om ee oe eee) 4 ~ Be ee in hh rr Pe tr) at hes © 80s Oe age 08 
oo ¥y abt coe os com oe . +s WhO eee ge, s i tee" oe SY oom og 8 sede 
ee ° es » “* om « © amOn. 6 re | ~ 8 eMVewe ae -m a 
trot) one Ste one Scie: ay " te ey - “e ww *8 30 Boeas, ay ao tw, 
oF Be ee we ve Bamenh eee. "© Se eo» ' » es = Buy ie . Lal en Se. mae e*e +0 
eh deh To > = © 0s de Seam « oe = e . . wa o- se eee, os — 2 oon 2 pares 
5 7 4 z . ee ee ars eee, oe Sey Gee ° Sy ek ete Les ’ ooh Pos His ah, we sbrone Ad ban] “Sw aioe wea e. 
Pose 5 oy © WY ee 6 mse. og © mace oe in) ae sate oP 40 noon’ ow 9 $0 % gd ant sof: 2 018s way one ‘ee  % aeuapen a eee: : 
wer moe . @tom a4 . magi Sie ere V ohee meting Ree Me eh eiccnces eee . =, an — Att ee wee 
ons "OR . ot = @ ae etpaces 6, ore 8 and one om Aen wanes g, = Shree net ge SO 08P0e Rae & 2 maw Cine be dendiie. Pete - 
OF Sar any Pe@eer os Fae s0ebe « eo teres “we , « 08s = im pe, eee ne ge e @, 6% © eases w@ wit «ls @o“eee Peet hs © @cs aceng oe Rene 
© beds Saree jew e so 8 dd ey ee we, Meee nen 2 > ote SS a i ematnigs Nomen ge meas verre wee ok . og ae NOt ten, Pye edhe tind 
POP 0800 Stigee * oie . ‘ ; a Oe ORAS ere hginetedere hee ne) Fe ey 7) 
erery "a OUSEe ws dtsee * et, sone we Be 9 = ome Oe Fe ites , oor ewe eras “er 13 e2e Gta © ee, Oo “Hayne ‘@u®. wid mew Oar. me wag pen Pantene. 
S gtte © ae ae hy are St cyte op al Tag anne » ome 1S 60 Wa ** osu, sone * SBatOteR Og oe ts y ata je) a Rag nes = Be wate & Chats cum we oean 
=~ Me ° 1 7e ot 96 ° § oye 9 ae ended «wp Coen 5 ve = Sen = OPeages. E WSe os geuy Wil se 002s cafe ~— Boome een) od newer ae Sr apy 
stee 8 more * be é Mew We [eth a tine soul wd so 8 MT te AN haley VD Sree alan eeth 8-80 “2g “mn eee wee a 0emt ee Gar ene 
fe Sry . hh * ‘ oo 2 . S ow webw wee BED MPCE 2D orag te 8 * ot %.” be iee Mme eae oe we 
aoe oom | w a2 : = <2 eh Oey arr, S\een awe Og anes 
Put ws * mT) e ah OO a © Se N we ot 65S o St Se ome 8% Ute es ag fee ethene ene eM ee WEP: Come ty ~an at 
eo Sewic sf ¢ 4 * @e % Bed athe r 7m 6 8 lg tm On rm o © @ a Ben asgns Hoe 5 ce bibs te 8 8) Se a seerne tO Meee wee 
° e * se. & . Sr ra, Woe ode my F oe Silitenel ccs iehere e Been sige VaPee.g He owe att 8 4 rs paid Ao we lS antes 
“se wee i tronee @ me. a a ont @ eng. © Owes wee taht TM FRMeme.  wntnte tei} ayer qrsd age Seen bed 
> we . & *g wren Sw © ewleem , * * Sscienuilies 6 ou be. i iw 8 oe 0 eee: ot" bp tae 5 ean unwes, 
fet, . fr ese op aie z “ =e et SSE ESSE BES ee 
é oe eu “ Ob ay ele Me- cede A rt er Parr © Pee e otk Risa.) & ras » ee tery Cate she Wwe Ge ao a Sse 
acta nea II odd nee dood cog bd a oe 08 Semmes se 8g ansy . com te bee 8, . Peres Sre, 00 © Ott, pee tie mY ne . " 7 ee nee 
oN etHteSs ee eecirlatecs Ss . 00 ‘eimce « ve Eo ae Heras c i . < : « a vte #0 Se Reb eo cart! r o ten eel twie ae gate me EO tts ~~ Sh" GR eae, 
Pomeny St mt no pre ead ° * = “Tyee 82 peste. 6 pad ee : * Bomn a cages Sot tx tigsyre os we hadon~gnan tee. 
vO 898 Q be cms bad le. J ° © * 2f Feee de gee ole . Me one © ape = On 0 8 ° eo S22 woe oe eR oe = Gee 
hs heP-oos 4. ny Cr ier @* Gee oe es om nee ° Peewee y © 0 ‘cee Ah 2 save 26" s : wk we pepe ege - ° “en. 
ROS ee Pee 2 cage cue as ees ee 4 1% sheet Sree SP sere aes = 
10 ts aw 6 Pewee, see e 48 *B0, = © ” of. = 5 hav se8. nd “ hit alnad aenntye 
© Ati. a mene a ABO ac: ek es ce fe tyeae te A onl. eee Heat ee Sy oe teow pate MY ABSA ES ES A tgeadag oe iss ee Se ee = 
5 @  wWe.boy .. Pome ay * Nw a ae “ect an . os oe 8 oy ee ouse mee ow. e Caer Ore me be oe ee 66 © Yetnee 2% ate we bs ote On SRG Oh 9 O50? omen af mg pede) | § tS we Cony Seed Reyes \ogms as O- egte> 
Sw loNin Ste Baye POOH Ee bee up Secours 6 re e tee Were 5 we euue wf ° « . 20% Toe we eat ett alemapacte So's 4 FR OG Repwts 8 Oe oe, Wirth lowe i oa, by emte-2 2 a LS PA patent 
10 e ° o = oar ®®.90.0 7} iS Som ee tee, ™ em? SO Pe Made. MO) we capers 28 .e  OPo * cures a= Sani 
1@e oe 2 5 Oh 2! fan 88 fom ee @5 7s em S we wn Aerts. ef % See see, ~ a ae Ao) ate ate POR Wy ne Sate 
were etm 200 a Bees 5 es ey SAMA on ng ene =o? 04s, 6 eh Ses %o tes e, & ons | CY oe | eee a 28 eee Gree tase *e0e te oes, oe § 2% tet cam Me wet Haig 50000 ogee’ or) a oe ps 
Weow oyee £8 Pw ee 8 * vane 8 + 8 av, . Sess e0 #8 \eer) c e ¢ » On ry Pte @_ogme. a8 - : ee As sera Pie e wenn tas * ‘ade sees OM 0.0 sede HS METAR <a eat SA Meh vbr wg re 
wr=e Neen Ulam, |g Heer og Se ab eas teh es ges x 5. ae ee ~~ se « wager. *8,! She at Tina tart. Px <2s eat.am coche P Oxo ig wren Sot he etree ey erypr eran, Wag en 
oe an ah Ch . © tet, * 4 “awe ~ ¥ Rete eee wees + paves Se 0 °e-.2 Oo. Ree cee wd Tees 6 aes ere Fe 
«Masher. reel Eee Oe or abe : bia AL CTE: nn ie he —_ z y oe Y Pah ae Se Let ty a oun. “we SmWM Somer g 5s coe MWA E SB opt ee S235 e ORE = 
. 6 eae =e ens © Lal) Ses ee ? Ad Miwpome ® es e ee «4, La « = ° een 
SOMA Re on ont,, © set? ooo < ft fh od We : . ote a . os 7 4 mad ts 2 Ore 5. de rn ee tos B+ e® ast. ag, Root te ke ace CAE, ‘sola Pee 4 eae om caspe 
ee ee ea the . ee eo wee Peete we ete ee . we = wee v PS omW wierslee «5 aus tet 84 Nee ‘ote oo enaT ern ew 0 win eons ae an Pet oe we: a SA ee aneers bad 
8 58-8" ete oy oe Senet Te we eS Po ows ae ail ® we ~ eu = S MAP ER ee, wg cane tt OoBee Water m6 tet S.vi wos Penta Sk Sc) AS Eee ‘ett Pek ~< Raat. 
° See Peay © 9O'"bem eree & Osh eoow ° se me tp ad SP ehete oe Soule ee ates com ote. ge ar] see reverdeat — 
".. s i Oe ee ors we, . CCMES. wy, ¢ © (« « Some 5 “mais Pm eX, “seals ee, ’ ere tm ff OO bom Que * CET ee fen mide et) sis: © A Ebony 
" > . feree,  . te Sa ee ow Yen 6 fem =a we vee deity ot salfe " « Bim om / WA cme eA nop "2am: oretr Pe. ae 9s, Apatite Teesrenas Ta % 29 eon we ees 
| : Perm Phe as vmee te bee wm tp arse a - & * Ts a BR ocd go a meres. 6 ed wy aie end Obi a +. en te eee oe Zee. ‘ 
i © aoe B, © cat thw GON te ose wm ms : ° mn = ea Sa TRes oad Re MBereP ing wade: te ote st) ut ee rete 
BIH we a Lit . 4 ro S08) ay oe 8 8 ™ peo ‘ BA Rewer © wie * ato. mw" % 270ee tem BO Wp" Dove “peed ees tees” 
ae B-) oe met lanes mel et ae ree de 4 So ebV gery, Ty eens Be SOO see w& @ Bodns 
ebdind SE. Me a for = 944 0 stuar, ss wt SS be% der f > @ we eye, “eens ot avon 
ort On Fer aog, cee OORT Iss. S826 08 Ber as, ne % Wve able CL 30) | wporen ligase wae 
Sud 4 the 2 one sc” te tbete ey 8 ba, ae e Oe 0h os ota, 28m’ 0s gw Bg e ore es 
oe 2b es 186) | ies o¢ +B eye “em tack @ ee “ $ o be wae we “Oeecgng *8s wetehe < “3-@ sagas 
OOH 0. 8 oe . ea © #08 6, Fi ow & = te Age o°°8 pase & ht.% sm, 7 
% alge o8em sae b cee &e eon 6 on. Ae © ate Fey ea most wets . of , iia J 
at . hd "Pe 4 & te Bee se = 
Vwe «+ ab p oe ay . . 9 t we 
8 ot... ape oo * 908) eee ° ¢ w te 
© 00 88hee; ow Bune ae 7 s ar 
ah ea PT . a so fe, 
= =) bey a ek eve t ¢ Saad 
a Or oh, doe we 
O° BB true sneele.8 
we mee fy Paw €, 
= — -. a 





